Draft version January 6, 2012 

Preprint typeset using WT^K. style emulateapj v. 8/13/10 



o 
a 



O 

u 

o 
is 



o 



THE WIGGLEZ DARK ENERGY SURVEY: GALAXY EVOLUTION AT 0.25 <Z< 0.75 USING THE SECOND 

RED-SEQUENCE CLUSTER SURVEY (RCS-2) 

I.H. Li 1 ' 2 '*, H.K.C. Yee 2 , Chris Blake 1 , Sarah Brough 3 , Matthew Colless 3 , Carlos Contreras 1 , Warrick J. 
Couch 1 , Scott M. Croom 4 , Tamara Davis 5 ' 6 , Michael J. Drinkwater 5 , Karl Forster 7 , David G. Gilbank 8 , 
M. G. Gladders 9 , Bau-ching Hsieh 10 , Ben Jelliffe 4 , Russell J. Jurek 11 , Karl Glazebrook 1 , Barry Madore 12 , 
D. Christopher Martin 7 , Kevin Pimbblet 13 , Gregory B. Poole 1 , Michael Pracy 4 , Rob Sharp 3 , Emily Wisnioski 1 , 

David Woods 14 and Ted Wyder 7 

Draft version January 6, 2012 

ABSTRACT 

We study the evolution of galaxy populations around the spectroscopic WiggleZ sample of star- 
forming galaxies at 0.25 < z < 0.75 using the photometric catalog from the Second Red-Sequence 
Cluster Survey (RCS2). We probe the optical photometric properties of the net excess neighbor 
galaxies. The key concept is that the marker galaxies and their neighbors are located at the same 
redshift, providing a sample of galaxies representing a complete census of galaxies in the neighborhood 
of star- forming galaxies. The results are compared with those using the RCS WiggleZ Spare-Fibre 
(RCS-WSF) sample as markers, representing galaxies in cluster environments at 0.25 < z < 0.45. 
By analyzing the stacked color-color properties of the WiggleZ neighbor galaxies, we find that their 
optical colors are not a strong function of indicators of star- forming activities such as EWQOII]) 
or GALEX NUV luminoisty of the markers. The galaxies around the WiggleZ markers exhibit a 
bimodal distribution on the color-magnitude diagram, with most of them located in the blue cloud. 
The optical galaxy luminosity functions (GLF) of the blue neighbor galaxies have a faint-end slope 
a of ~ —1.3, similar to that for galaxies in cluster environments drawn from the RCS-WSF sample. 
The faint-end slope of the GLF for the red neighbors, however, is ~ —0.4, significantly shallower than 
the ~ — 0.7 found for those in cluster environments. This suggests that the build-up of the faint-end 
of the red sequence in cluster environments is in a significantly more advanced stage than that in the 
star- forming and lower galaxy density WiggleZ neighborhoods. We find that the red galaxy fraction 
(fred) around the star-forming WiggleZ galaxies has similar values from z ~ 0.3 to z ~ 0.6 with 
fred °° 0.28, but drops to f re d ~ 0.20 at z^0.7. This change of f re d with redshift suggests that there 
is either a higher rate of star-forming galaxies entering the luminosity-limited sample at z^0.7, or a 
decrease in the quenching rate of star formation at that redshift. Comparing to that in dense cluster 
environment, the f re d of the WiggleZ neighbors is both considerably smaller and has a more moderate 
change with redshift, pointing to the stronger and more prevalent environmental influences on galaxy 
evolution in high-density regions. 

Subject headings: galaxies: evolution — galaxies: photometry — galaxies: luminosity function, mass 
function 
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1. INTRODUCTION 

The subject of galaxy evolution has been widely stud- 
ied using both photometric and spectroscopic data over 
wide redshift ranges. In the Northern sky, the SDSS sur- 
vey has mapped out a vast region of the nearby Universe, 
and numerous studies have investigated galaxy proper- 
ties and environmental influe nces using its da ta. Other 
survey s, such as COMBO- 1 7 flWolf et qMpI) and COS- 
MOS (jScoville et al.l 120071 : iLillv et al.ll2007D . have spent 
much effort to explore galaxy properties and evolution 
in the more distant universe, out to redshift ~1 and be- 
yond. It has become clear that galaxy colors exhibit 
a bimodal distribution at all redshifts to at least z ~ 1, 
with a relatively narrow red sequence dominated by non- 
star-forming gala xies and a blue cloud of star-forming 
galaxies (e.g., |Strateva et al.l 120011; iBlanton et al.l [2003; 
iBell et all 120041 : iWillmer et al.l I2006D . The fraction of 
red-sequence galaxies (or blue cloud galaxies) changes in 
different environments and at different redshifts. It has 
been found that red passive galaxies tend to populate 
dense environments and blue star-formin g galaxies are 
more common in less dense regions (e .g., iDressler et al.l 
119801 : ICooper et al.lEHMlLi et al.ll2009D . In clusters from 
z < 0.1 to z ~0.5, the fraction of blue galaxies in- 
creases from a few per cent to ~30%, and reaches ^70% 
at z ~1 (e.g., iButcherlF Oemler 198 j; lLoh et al.l 1 2008: 
iMahaian fc Ravchaudhur v 2009; Haines et all 12009( 1. 

It is believed that the red sequence in galaxy 
clusters is assembled from the top down, being al- 
ready largely in place at the bright end by z ~ 
1, with the fa i nt-end filled in at a later time (e.g. , 
Bell et al. 2004; Tanaka et al. 2005; Willmer et al. 2006; 
Stott et all 120071 : iGilbank et all 120081 : iDe Lucia etaLl 
2009h . Since there are relatively fewer stars formed 
in red-sequence galaxies, the build-up of the red se- 
quence has been argued to be driven by the global sup- 
pression of star formation through environmental-related 
processes, such as galaxy merging, galaxy harassment, 
gas s tripping, or gas co nsum ption by star-forming disks 
(e.g.. IDressler fc Gunnl 119831: iBarnes fc Hernquistl [l99lt 
iMoore et al.l ll996f T The build-up of the red sequences 
in galaxy clusters can be considered as the gradual 
loss of late-type prog enitors over a Hubble time (e.g., 
Ivan Dokkum fc FraiuJl200H : iKavirai et al.ll2005f ). 

In field environments, galaxies also exhibit a bimodal- 
ity in their color distributions, and form a red sequence. 
However, the fraction of field galaxies on the red sequence 
is greatly lower than that in clusters. Nevertheless, 
there is also a deficit of faint red-sequenc e field galax- 
ies, both at higher redshift (e.g ., at z ~0.8. iTanaka et all 
I2005L IgOOgUW einer et al. 2 00 j) and at the present day 
(e.g.. IWvder et al.l 120071 : lBlantonll2006D . indicating that 
the assembly of the red sequence is still incomplete in 
low-density environments. As the majority of star for- 
mation at all redshifts is contributed by blue late-type 
galaxies, it is interesting to probe galaxy evolution from 
the perspective of blue star- forming galaxies. Especially 
the average star-formation density is evolving rapidly 
with redsh ift in the field at least a factor of ~ 10 since 
z ~1 fe.fr. iMadau et ai1ll998t iHopkins fc Beacoml 120061 : 
IGilbank et aTll20lol ). 

Many studies of galaxy evolution beyond the local 
Universe focus on red galaxies or the cluster envi- 



ronment (e.g. iBalogh et al.l 119971; iGladders et al.l 119981 : 
iKodama et aLll!998b iLemaux et al1l201Q[ ). While the ev- 
idence is clearer on the question of cluster environmental 
influences, the star formation of galaxies in field regions 
is still ambiguous. Because star formation is still active 
in blue galaxies, they provide a more direct observation 
on the actual dependence of star formation rate on envi- 
ronment and also its evolution. 

In this paper we will use a combination o f the Sec- 
ond Red-Sequence Cluste r Survey (RCS2; (jYee et al.l 
120071 : IGilbank et all 1201 1ft) and the WiggleZ spectro- 
scopic survey (jDrinkwater et al.l 12010) to study galaxy 
evolution up to z ~ 0.7. This combination pro- 
duces one of the largest photometric and spectroscopic 
databases at intermediate redshifts, covering a total of 
~300 square degrees with g'r' z' photometry and optical 
spectra (4700A-9500A) for -120,000 blue star-forming 
galaxies at 0.2 ^z^l. Such a combined data set provides 
a great opportunity to explore properties of the galaxy 
population and its evolution at the intermediate redshift. 
The WiggleZ spectroscopic survey targets primarily blue 
star- forming galaxies, using UV fluxes as its main selec- 
tion criteria, along with a set of complex optical selection 
rules (see $2]). While the WiggleZ spectroscopic catalog 
provides a large sample of star-forming galaxies cover- 
ing a significant redshift range, its complex optical selec- 
tion criteria make its direct application for investigating 
the evolution of star-forming galaxies complicated, if not 
impossible. However, they provide a valuable database 
as a catalog of markers of regions where star formation 
is prevalent, which are likely low galaxy d ensity regions 
of the universe. Inspired by the work of lYee fc Greenl 
(|1987| ). who used low-redshift quasars as markers to de- 
rive statistically the luminosity function of galaxies asso- 
ciated with quasars, we approach the task by exploring 
photometric properties of the galaxies around WiggleZ 
galaxies, which provide an unbiased census of galaxies in 
regions of strong star formation. The RCS2 survey pro- 
vides complete and relatively deep optical photometric 
catalogs for half of the fields used by the WiggleZ survey, 
and they are used for the analyses of the WiggleZ galaxy 
neighbors. Our work shows that probing the properties 
of the neighbor galaxies of markers statistically can offer 
a powerful method in studying galaxy evolution. 

The structure of the paper is as follows. We describe 
briefly the WiggleZ spectroscopic project and the photo- 
metric RCS2 survey in $21 Our method of constructing 
color-color-magnitude cubes of the galaxies associated 
with the WiggleZ galaxies is detailed in $3j The results 
are presented in $4j where we probe the color-color plots, 
color-magnitude diagrams, luminosity function, and red- 
galaxy fractions of the neighbor galaxies. We discuss the 
results in $5] and summarize our work in $6] We adopt a 
cosmology of f2 m =0.3, ^a=0.7, and iJo=70km/s/Mpc. 

2. THE SURVEYS AND DATA 

The basic assumption used in this work is that, because 
galaxies cluster, excess galaxies counted around a marker 
of known redshift are in the same redshift space as the 
marker, allowing us to measure their intrinsic photomet- 
ric properties such as luminosity and rest-frame colors. 
To this end, we require a sample of galaxies with spec- 
troscopic redshifts and photometric data of the complete 
field of the spectroscopic sample. In this section, we first 
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describe briefly the WiggleZ and RCS2 surveys, which 
provide the spectroscopic and photometric data, respec- 
tively; we then present the actual sample of the WiggleZ 
markers used. We also present a comparison sample of 
markers, obtained as part of the WiggleZ observing runs, 
based on positions of RCS2 clusters. 

2.1. The WiggleZ Spectroscopic Survey 

2.1.1. Target Selection 

The WiggleZ Dark Energy Survey is a spectroscopic 
survey of 240,000 UV-selected emission-line galaxies, de- 
signed to map a cosmic volume of ~lGpc 3 . Its primary 
goal is to precisely measure the scale of baryon acous- 
tic oscillation (BAO) imprinted on the spat ial distribu- 
tion of these galaxies at 0.2 < z < 1 (e.g., IB lake et"aTl 
20111). The details o f the survey are presented in 
Drinkwa ter et al.l (|2010f ); here, we present a brief sum- 
mary. 

The survey selects targets from areas totaling ~ 1,000 
deg 2 from 7 equatorial regions. Target galaxies are se- 
lected using the FUV and NUV data from the GA LEX 
Medium Imaging Survey (MIS; iMartin et all 120051 ) us- 
ing the criteria of FUV - NUV > 1 or no FUV detec- 
tion. The targets must also satisfy NUV < 22.8 and the 
NUV signal-to-noise ratio S/N > 3. Further selection 
criteria based on optical photometry are also applied to 
attempt to maximize the probability that the targets are 
blue star- forming galaxies at z > 0.5, the primary sam- 
ple for the WiggleZ project science. The optical data 
are o btained from SDSS DR4 (Ad elman-McCarthy et all 
I2006D and RCS2 (see g22]). In order to select blue star- 
forming galaxies and exclude spurious matches between 
GALEX and optical data, all WiggleZ targets must also 
have —0.5 < NUV — r' < 2. To avoid low-z galax- 
ies, a 20 < r' < 22.5 criterion together with two dif- 
ferent sets of optical color-color selections are applied. 
All these selection criteria give a target density of ~350 
galaxies/deg 2 , or, ~2.6±0.2% of optically detected galax- 
ies. There is no further morphology selection to remove 
any 'stellar' objects, since galactic stars ought to fail the 
survey selection criteria. 

The SDSS data have a depth of [22.0, 22.2, 22.2, 21.3, 
20.5] in the u'g'r'i'z' passbands, while the RCS2 data 
have a much deeper depth, with average 5a point source 
limits of [24.4, 24.3, 22.8] in the g'r'z' passbands. Since 
we want to study the galaxy population properties and 
their evolution to as high a redshift as possible, in this 
paper we use only the RCS2 regions of the WiggleZ sur- 
vey. 

2.1.2. Observation and Data 

The WiggleZ observations were conducted using the 
AAOmega spectrograph (the former 2dF upgraded; 
iSharpI 120061 ) on the 3.9m Anglo- Australian Telescope 
(AAT) from Aug 2006 to Jan 2011. AAOmega is a 
fiber spectrograph containing 400 fibers including 8 guide 
fibers. Each fiber has a diameter of 2" . The field of view 
is 2 degrees in diameter. The typical exposure time is 
60 min per AAOmega configuration. This exposure time 
is too short for allowing a significant detection of galaxy 
continuum for the fainter galaxies, but sufficient to de- 
tect emission lines for redshift measurement. Using the 
580V and 385R gratings for the blue and red arms with 



the 670nm dichroic, the spectra have a wavelength range 
from 4700A to 9500 A, with a dispersion of ~l.lA/pix 
in the blue arm and ~1.6A/pix in the red arm, providing 
spectral resolutions of ~ 3.5A and ~ 5.3A, respectively. 
The observing conditions varied significantly, with the 
seeing typically ranging from 1-2. 5". 

Detailed descriptions of the data reduction tech nique 
and reliability are given in Drink water et al.l ([2010) and 
summarized here. The data were reduced during each 
observing run using the automated 2dFdr software de- 
veloped at the Australian Astronomical Observatory, in- 
cluding bias subtraction, flat field, and wavelength cal- 
ibration. The redshift of each spectrum was measured 
using an evolved versio n of runz, which w as the soft- 
ware used for 2dFG RS (jColless et aLll2QQlh and 2SLAQ 
(jCannon et al.ll2006l ). The software has been modified to 
optimize the use of emission lines to derive redshifts. The 
commonly detected emission lines in WiggleZ spectra are 
[OII]A3727, Hp, [OIII]A4959/5007, H a , and [NII]A6583. 
Even though runz automatically generates an integer 
quality flag (Q ZS pec) in the range of 1-5 based on how 
well the template fits a given spectrum, all WiggleZ spec- 
tra were extensively checked visually, and each spectrum 
was manually assigned a new quality flag. The redshift 
confidence increases with larger Q zspec . The redshift re- 
liability has been cross-checked internally using repeated 
galaxies. A subset of redshifts were also compared to 
DEEP2 galaxies. While there may be some debates in 
distinguishing Q zspe c =4 and Q zspe c=5 objects, as they 
are close to being 100% reliable, the critical separation 
occurs between Q zspe c=^ and Q zsp ec = 3. For objects 
with Q zsp ec=3, the redshift reliability is ~79%. 

As of Oct 2009, the survey yielded -260,000 spectra 
in total from all 7 equatorial regions, and —160,000 of 
them are useful with Qzspec=3,4,5. Part of the Wig- 
gleZ spectral database has been released to the public 
at [http://wigglez.swi n.edu.au/dsl M ore details about 
the pr oject and data can be found in iDrinkwater et all 
(|2010h . 

2.2. RCS2 Photometric and Random Catalogs 

The RCS2 is a — 1000-square-degree imaging survey 
in z\ r' , and g f with the goal of identifying a large 
sample (> 10 4 ) of galaxy clusters up to z — 1 for the 
purpose of constraining cosmological parameters using 
the galaxy cluster mass function and studying galaxy 
evolution. The survey was carried out using the one- 
square-degree camera MegaCam at CFHT. The survey 
targets 12 regions of sky with areas varying between 36 
to 100 square degrees. About half of the WiggleZ fields 
use RCS2 positions and photometry for target selection. 
The three-color photometric catalogs of galaxies in these 
targeted fields are used for our analysis of companions 
associated with WiggleZ galaxies. The details for the 
survey and photometric catalog production are described 
in Gilbank et al. (2011); here, we provide a very brief 
summary. 

The RCS2 photometric catalogs are created using an 
automated pipeline, with algorithms for object finding, 
photometry, and star-galaxy classification based on those 
from the program Picture Processing Program (PPP) of 
lYeel ([19911 ). The high -precision photometric catalogs in 
g'r' z' are calibrated using the colors of the stellar locus 
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Fig. 1. — Zspec distribution of 62785 WiggleZ galaxies with red- 
shift quality flag Q ZS pec=3,4,5 at z < 1.5 from four RCS2 patches. 
The redshift binsize is Az=0.01. The WiggleZ selection function 
is designed to optimize observations of high-redshift galaxies thus 
a significant number of galaxies with z ~0.3-0.6 are removed from 
the sample. In this paper we focus on the redshift range between 
z=0.25 and z=0.75, which are marked by the vertical dashed lines. 

combined with overlapping Two-Micron-All-Sky Survey 
(2MASS) photometry. This technique yields an abso- 
lute accuracy of better than ~ 0.03 mag in colors, and 
~ 0.05 mag in the r'-band, verified via regions that over- 
lap with the SDSS. The survey reaches average 5<r point 
source limiting magnitudes for z',r' ,g' of 22.8, 24.3, and 
24.4, respectively, approximately 2 magnitudes deeper 
than the SDSS. Absolute astrometric calibration is accu- 
rate to better than 0.3". 

A key feature in using the RCS2 catalogs is the avail- 
ability of random catalogs. These are random points gen- 
erated to populate the survey area of each region with a 
uniform density of one per 10 square arcsecond. These 
random points allow one to map out the areas where 
there are no data, including chip gaps, bad columns, 
bright star halos, saturated pixels, meteor trails, and 
other cosmetic defects. These catalogs are crucial in es- 
timating the area sampled by the data, as in generating 
proper background count estimates in the analysis per- 
formed in this paper. 

2.3. The Sample of the Markers 

Since the RCS2 imaging is much deeper than the 
SDSS data, the marker galaxies are chosen from the four 
WiggleZ-RCS2 regions. They are RCS2 0047+00, 0310- 
14, 2143-00, and 2338-09; namely, the Olhr, 03hr, 22hr, 
OOhr fields in the WiggleZ survey layout. The use of the 
RCS2 regions allows us to measure photometric proper- 
ties to a much higher redshift. 

We use the WiggleZ data taken prior to Oct 2009. 
There are 62785 spectra in total with redshift quality 
flag Q 

zspec ^ 3 from these four WiggleZ-RCS2 regions. 
The redshift distribution is presented in Fig. [TJ The me- 
dian redshift is z ~0.59. We focus on the redshift range 
between 0.25 and 0.75, giving a total of 41041 spectra. 
The lower z = 0.25 boundary is chosen so that our bluest 
passband g' is still blueward of the 4000A break at the 
lowest redshift bin, while the upper redshift limit is set 
based on having sufficient depth in the imaging data for 
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Fig. 2. — The observed g'-r' colors as a function of redshift for 
WiggleZ galaxies. Only 1/5 of the sample is plotted for clarity. 
The red, green, and blue solid curves are the colors from GISSEL 
evolving spetra described in §2.3, while the corresponding (dot- 
) dashed curves are the colors for the non-evolving spectra used 
for k-correction. A 'hollow' feature with a deficit of galaxies is 
observed at z=~0.3-0.6 and g'-r' ~0.7 due to the 'low-z rejection' 
in the survey selection criteria. The red dots are galaxies that 
satisfy the 'low-z rejection' criteria but were still observed during 
early stages of the WiggleZ survey (before these criteria were fully 
implemented). 

the analysis of the companion galaxies. 

The observed g'-r' colors of the WiggleZ galaxies as 
a function of redshift are presented in Figure El We 
also o yer plot three dust- free mo dels generated from GIS- 
SEL (jBruzual fc Charlotl I2QQ3I ) . The red dashed curve 
is the S ingle Stellar Popu lation (SSP) model using the 
Padova (jBertelli et al.|[l994[ ) e volutionary trac ks with so- 
lar metallicity Z=0.02 and the lChabrierl (j2QQ3h IMF with 
a zero-redshift age of 13 Gyrs, representing evolved early- 
type red galaxies. The green dot-dashed curve is a r 
model with an exponentially decreasing star-formation 
rate with r = 1 Gyr, generated with a metallicity of 
Z=0.0001 and a zero-redshift age of 11 Gyr, represent- 
ing mildly star- forming spiral galaxies. The cyan short- 
dashed line is a model made with a constant star forma- 
tion rate, representing star- forming late- type/irregular 
galaxies. 

We observe two features from Figure [2] First, red ellip- 
tical galaxies are absent from the WiggleZ sample. Most 
of the galaxies populate the region between the r- and 
constant star-formation models. This is expected, be- 
cause the WiggleZ galaxies are selected by the UV flux 
to be star-forming galaxies. 

Second, there is a 'hollow' region with a deficit of galax- 
ies at g'-r' ~0.7 between z ~ 0.3 — 0.6. This is also mani- 
fested as a dip in the redshift distribution of the WiggleZ 
markers at 0.3 < z < 0.6 in Fig. [TJ This 'hollow' feature 
arises artificially due to the 'low-redshift rejection' (LRR) 
criteria based on g' — r' > 0.6 and r' — z' < 0.7 (g' — r') in 
the WiggleZ 's selection criteria used in the RCS2 regions, 
in the attempt to maximize z > 0.5 galaxies. Thus, the 
LRR criteria actually remove galaxies at redshift up to 
z ~ 0.6, producing the significant broad dip in the red- 
shift distribution of the WiggleZ galaxies. Nevertheless, 
the 'hollow' region has more galaxies in it than expected 
with the LRR criteria. There are 8597 galaxies in the 
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sample that actually meet the LRR criteria but are still 
included in the WiggleZ sample. This is because the LRR 
criteria were developed after the early observing runs and 
thus some galaxies which were initially observed would 
have been rejected later on by the refined selection crite- 
ria. We find that these galaxies, overplotted in Fig. [2] as 
red dots, are primarily distributed over 0.25 < z < 0.65. 
Thus, while the WiggleZ survey intends to target blue 
star- forming galaxies, an examination of Figure [2] indi- 
cates that the sample is composed of a range of star- 
forming galaxies, with colors consistent with starburst 
to constant and mildly star- forming galaxies. 

Finally, we note that any direct comparison of the 
properties in this spectroscopic sample as a function of 
redshift is not straightforward, due to the 20 < r' < 22.5 
selection criterion. This criterion produces galaxy sam- 
ples of different absolute magnitude ranges at different 
redshifts. 

2.4. The RCS2 Cluster WiggleZ Spare-Fibre 
(RCS-WSF) Sample 

During the WiggleZ observing runs, a small number 
of AAOmega fibres were used for targets from different 
projects, using identical observation parameters and data 
reduction techniques. In addition to the sample of the 
UV-selected WiggleZ galaxies, there are ~3000 spectra 
targeting RCS2 cluster galaxies as part of the WiggleZ- 
RCS2 collaboration. These galaxies are selected from a 
preliminary sample of RCS2 clusters at z < 0.5, chosen 
as possible high ranking, bright, red-sequence galaxies in 
the clusters. Most of these galaxies are at z < 0.45 and 
have a mean redshift of z ~0.28. The details and scien- 
tific results using this WiggleZ-RCS2 cluster subsample 
will be presented in a future paper. For the purpose of 
this work, they serve as an excellent comparison sam- 
ple of markers to the WiggleZ galaxies, as they are red 
galaxies in dense environments. Since these galaxies do 
not cover the same redshift range as the WiggleZ galax- 
ies, the comparison is only available at lower redshifts. 
We will refer to this sample as the RCS2 WiggleZ Spare- 
Fibre, or RCS-WSF, sample for the remainder of the 
paper. 

3. METHOD 

With the assumption that the WiggleZ marker galaxies 
and their neighbors reside in the same spatial regions, we 
can construct a net color-color-magnitude (CCM) cube 
of the neighbor galaxies, so the colors and magnitude 
information of the neighbor galaxies can be preserved. 
The xyz axes of the cube represent r' — z' , g' — r', and 
r', respectively. Thus, the net counts as a function of 
luminosity in the r' passband, for instance, can be com- 
puted by summing the values in the x and y axes along 
the z axis. Essentia lly, we adopt the me thod used in 
IGilbank et all ([2008) and lLoh et al.l IpQOl for creating 
color- magnitude diagrams, but extend the concept to a 
3D cube. To produce the net CCM cube, we subtract a 
background CCM cube from the total-count CCM cube. 
The cubes are made in both observed and rest frames. 
We detail the methods below. 

3.1. Observed Color- Color- Magnitude Cubes 

To construct a CCM cube, we first identify all galaxies 
in the RCS2 photometric catalogs with r' < 24.0 within 



a projected comoving radius of r p Mpc to a WiggleZ 
galaxy. The WiggleZ galaxy itself is excluded in this 
process. These galaxies are namely the 'neighbors' to 
the WiggleZ galaxy, and their observed r 1 — z' \ g' — r 1 \ 
and r' are gridded into a cube with a binsize of 0.05, 0.05, 
and 0.1 along the xyz (color-color-magnitude) axes. The 
CCM cube of the control field (i.e., the background) is 
constructed using all galaxies with r' < 24 in the same 
RCS2 patch. A single RCS2 patch is sufficiently large 
(typically 81 square degrees) to provide excellent back- 
ground statistics, and by using the same patch, it also 
ensures a minimal systematic effect. The control cube 
has the same binsize as the cube of the marker neigh- 
bors. The count of each element in this control cube 
is then scaled by N ran ^ n /N ran jot, which is the ratio of 
the number of the random points within the aperture 
(N ranjin ) to the total count (N rarijtot ) in the patch (see 
The typical N ranM /N rariyto t is - 2 x lO" 6 . A net 
CCM cube is obtained by subtracting the scaled control 
cube from the cube of the neighbor galaxies, i.e., 

CUbe ne t — CUbe ne ighbor CU-bebackgroimd ^ -^ran,in/N ran? tot • 

These net CCM cubes from individual markers can then 
be stacked to form the total CCM cube. 

3.2. Rest-Frame Color- Color- Magnitude Cubes 

One approach to obtain a rest-frame CCM cube is to 
convert it from an observed one which has been described 
above. However, it requires a large amount of computing 
time to k-correct each element of an observed CCM cube. 
An alternative is to compute the rest-frame magnitude 
and colors of each galaxy first, then construct the rest- 
frame CCM cube using the same procedure as building 
the observed cube. Since there is no actual redshift infor- 
mation for the neighbors, it is assumed that all surround- 
ing galaxies are at the same redshift as the marker. A 
control field cube is computed for each marker, for which 
we also convert all galaxies into "rest-frame photometry" 
using the redshift of the WiggleZ galaxy. 

The k-correction is derived using tables gene rated for 
each of the GISSEL (jBruzual fc Chariot! l2003h models 
described in §2.31 Each table contains galaxy colors of 
the model and the k-correction values for each passband 
as a function of redshift. For a galaxy at a fixed red- 
shift, we derive the k-correction using the model grid 
by interpolating (or extrapolating in some instances) the 
model colors to match the observed galaxy colors. The 
model colors here are the observer-frame colors of GIS- 
SEL galaxies with non-evolving spectra, which are over- 
plotted as curves in Figure O We use the g' — r' color 
to derive k-corrections for the g' and r' passbands, and 
r' — z' for the z' magnitude. We have compared our 
k-correction results to the SDSS galaxies in one region, 
where their k-corrections are available from the official 
SDSS database. Our method in deriving the k-correction 
yields a good correlation with the SDSS values. We 
use z = 0.5 as our reference redshift, and all the rest- 
frame photometry is computed relative to this redshift 
as 5 M = m - DM — (k — £^=0.50), where DM is the 
distance modulus. For reference, Figure [3] plots the k- 
correction as a function of redshift for each filter for the 
three spectral types. 

4. RESULTS 
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Fig. 3. — k-correction as a function of redshift for g'r'z' niters. 
The dashed red, dot-dashed green, and short-dashed purple curves 
are generated using the same model in Fig. [2] which represent 
evolved early-type red, mildly star-forming spiral, and star-forming 
late- type/irregular galaxies respectively. 



TABLE 1 

Galaxy Counts 



redshift 


N WiqqleZ 


N net 


■^backqround 


0.25-0.35 


6885 


10332.0 


118828.0 


0.35-0.45 


6090 


7811.85 


61085.1 


0.45-0.55 


6857 


8021.05 


46571.0 


0.55-0.65 


10869 


9948.05 


54039.0 


0.65-0.75 


10340 


9084.79 


40491.2 



a Note: 
radius r 



N ne t and Nb ac k groU nd are within a projected co-moving 
p =0.25Mpc. 



In this section we derive the various photometric prop- 
erties of the galaxy population associated with the Wig- 
gleZ marker galaxies. Our sample is limited to 0.25 < 
z < 0.75, and is divided into five redshift bins with 
Az=0.10. We refer to these redshift bins as z = 0.3, 0.4, 
0.5, 0.6 and 0.7. We group the WiggleZ galaxies into the 
redshift bins, and stack all the net CCM cubes within 
each redshift bin. The numbers of the markers and their 
neighbors in each redshift bin are listed in Table [TJ A to- 
tal of 45,198 net galaxy counts around 41,041 markers are 
used in our analysis. All the CCM cubes are made using 
galaxies within a projected comoving radius r p =0.25Mpc 
from the markers. The choice of this r p will be justified 
in §4.2. The CCM cubes of the RCS-WSF sample pre- 
sented in §4.5 are made with an angular-diameter radius 
of r p =0.25Mpc instead of a comoving one, since cluster 
galaxies are considered gravitationally bound, although 
the results are similar when using either a comoving or 
angular-diameter radius due to their redshift range. The 
galaxy counts in the RCS-WSF sample are tabulated in 
Table [2j The average number of net companions to the 
markers in this sample are about an order of magnitude 
larger than that for the WiggleZ sample. 

4.1. Random Marker Fields 

To test the reliability of background subtraction in our 
method, we construct CCM cubes based on the positions 
of 4000 randomly drawn points from the random catalogs 
in each of the four RCS2 patch. These random points 



TABLE 2 

Galaxy Counts for the RCS-WSF Sample 



redshift 


Nrcs-wsf 


N net 


Nback ground 


0.25-0.35 


416 


8968.30 


12270.7 


0.35-0.45 


294 


4760.45 


5754.54 



a Note: 
radius r % 



Nnet and Nb ac kg roU nd 

=0.25Mpc. 



are within an angular-diameter 
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Fig. 4. — The N ne t of the random catalogs computed using 
r p =0.25Mpc. a): distributions of the N ne t for the 4000 random 
markers in each RCS2 patch, b), c), d): The mean N ne t per bin 
(Ar' = 0.1, A(g'-r') = 0.05, or A(r'-z') = 0.05) along each axis 
of the stacked cube over 0.25 < z < 0.75. No systematic offset is 
observed. 

are then assigned a redshift between z=0.2b and z=0.75. 
This gives us ~3200 random points in total for each red- 
shift bin. All the CCM cubes are built in the same way 
as described in ^3] but with the marker position and red- 
shift replaced. Because these markers are randomly cho- 
sen and not based on actual positions of any galaxy, we 
expect the average net excess in the neighbor counts to 
be zero when these net CCM cubes are stacked, if the 
background subtraction is properly handled. 

We compute a net neighbor count, N netl for each ran- 
dom point by summing the intensity in all elements of an 
observed-frame net CCM cube built with r p =0.25Mpc. 
The N net distribution together with its dependence on 
galaxy magnitudes and colors are plotted in Fig. 3J 
Panel a) shows the N net distributions of the 4000 ran- 
dom markers in each RCS2 patch. These distributions 
are statistically identical for the different patches. Sum- 
ming these distributions gives a median N net = —0.344 
and a mean N net = 0.032 ±0.027 galaxies, where the un- 
certainty is the rms of the mean. Thus, the mean of the 
net counts around random points is consistent with being 
zero. The relatively large negative value of the median 
of the net counts is the result of galaxies being clustered 
even on the projected sky, which results in a skewed his- 
togram of the net count distribution. Because there is 
no observed offset in < N net > among different redshift 
bins, we stack all the cubes over 0.25 < z < 0.75, and 
project the stacked cube along an axis of r', g' — r' , or 
r' — z' . The total N net of the stacked cube along each axis 
is presented in the b), c), and d) panels in the Figure. 
The < N net > counts are also not a function of mag- 
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Fig. 5. — The differential mean N ne t in each annulus without 
(top) and with (bottom) area normalization. In each panel, the 
solid red, dotted green, dashed blue, dot-dashed cyan, and dot- 
dot-dashed pink curves represent our redshift bins in an increasing 
order. There are ~1.5 net excess galaxies within r p =0.25Mpc, 
giving a surface density of ~6 gal/Mpc 2 . The surface density drops 
quickly to ~0 gal/Mpc 2 at r p > 1.5Mpc. 

nitude and colors, and have means of essentially 0, in- 
dicating that the background contamination is correctly 
subtracted, statistically speaking. Given these results, 
we are confident of our method in background correction 
and constructing the CCM cubes. 

4.2. Net Excess Galaxy Surface Density 

The WiggleZ survey targets blue star-forming galax- 
ies with a set of complex selection functions. Blue 
star-forming galaxies are believed to populate less dense 
environment compared to red passive galaxies (e.g., 
iDressler etl!lll980l : IWeinmaiin et aD2006[1CooDer et all 
l2007h . To investigate the characteristics of the neighbor- 
hood of WiggleZ galaxies, we probe the total net neigh- 
bor counts, N netl as a function of radius centered at each 
WiggleZ galaxy. The observed CCM cubes in a series 
of annuli are computed, and N net in an annulus is the 
sum of the intensity of all elements in the cube. Even 
though all the observed CCM cubes are constructed us- 
ing galaxies to a fixed apparent magnitude of r / =24.0, 
we note that the comparison among different annuli at 
a fixed redshift bin is still meaningful. Direct compar- 
isons among different redshift bins, however, cannot be 
made because the cubes are not limited to the same abso- 
lute magnitude depth for the different redshift bins. The 
mean N net in each annulus for the different redshift bins 
are plotted in Fig. [5] as a function of r p . The number 
of net excess galaxies within an annulus is not large, the 
maximum being only ~1.5 within r p =0.25 Mpc. Nor- 
malizing the N ne t by the aperture size, the mean surface 
density is a strong function of radius, being ~6 gal/Mpc 2 
within r v — 0.25Mpc and then decreasing rapidly with 
increasing radius and reaching ~0 at r p > 1.5 Mpc. Be- 
cause most N net excess is observed within 0.25 Mpc, we 
therefore use r p =0.25 Mpc to construct the observed- 
and rest-frame CCM cubes for our further analysis. 

4.3. The WiggleZ Galaxies and their Neighbors 
4.3.1. Observed Color-Color Diagrams 



Observationally speaking, galaxies appear to be pri- 
marily divided into two classes. One is red passive galax- 
ies and the other, blue star- forming galaxies. These two 
classes of galaxies form the so-called 'red sequence' and 
'blue cloud' in color-magnitude space. In fact, red galax- 
ies may be a mix of truly old passive galaxies and dusty 
star- forming galaxies, and they cannot be dist inguished 
well using a single optical color. However, IWolf et all 
(2005J) showed that dusty star-forming galaxies can be 
well separated from old passive ones in a color-color 
space, as long as one color brackets the 4000 A break and 
the other is at a longer wavelength. They found that 
dusty red galaxies actually form a continuous tail ex- 
tending from the blue cloud, while old red galaxies form 
a separate structure of their own (the red sequence). This 
makes the color-color diagram a powerful diagnostic tool. 

Figure [6] presents the observed color-color diagrams 
for neighbors in 5 redshift bins, where the colors of 
the three models from Figure [2] are overlaid as crosses 
for reference. For better visual presentation, the pix- 
els (A((/ — r') x A(r' — z')) in the color-color intensity 
plot are subdivided by a factor of 4 into smaller pixels in 
units of 0.0125 mag, and then smoothed by a kernel of 
10 x 10 pixels. The intensity scale is in units of counts per 
small pixels after normalizing the net counts to 1 x 10 4 in 
each redshift bin. We also overplot the WiggleZ galax- 
ies as the non-filled contours with a A(r'-z')=l offset 
for clarity. We observe that, for all redshift bins, both 
the WiggleZ galaxies and most of their neighbors pop- 
ulate a similar color-color plot, with the exception that 
the WiggleZ galaxies do not show a clump of red-passive 
galaxies. They both exhibit a continuous sequence in all 
redshift divisions. The sequence runs from the blue star- 
forming regions toward the red passive area, marked by 
the model colors; but no WiggleZ markers have colors as 
red as the red passive galaxies, indicating that they con- 
tain little dust. We note that this is likely a reflection of 
the survey design, as the WiggleZ galaxies are selected 
primarily by UV fluxes. 

Although most of the neighbor galaxies reside in the 
star- forming sequence, some neighbors populate the re- 
gion of passive red galaxies. These red neighbors are 
red-sequence galaxies, and we will discuss their proper- 
ties later in the paper. Some neighbors in the lower- 
redshift bins, however, exhibit redder (z f — r') colors than 
that expected from passive red galaxies at their fixed red- 
shift, and form roughly a continuous sequence from the 
blue star- forming galaxies. These neighbors are possible 
dusty galaxies. We note that such estimation is approx- 
imate, since the regions in the color-color diagram for 
dusty reddened star-forming galaxies may change at dif- 
ferent redshifts due to the shifting of the 4000A break 
in observed frame. At z ~ 0.6 and beyond, the use of 
the current color-color diagram to distinguish between 
dusty star-forming and red passive galaxies is not op- 
timal, because the 4000A break is shifted beyond the 
center of the r' passband, and having only one passband 
(V) at longer wavelengths is not sufficient to distinguish 
between passive and dusty SEDs. Using a rough color 
cutoffs to define dusty galaxies as (r r — z') — Ce > 0.3 
and g f — r f > AC where Ce is the color of the red el- 
liptical model in Figure [2] and AC is the color halfway 
between the red and green models in the same Figure, we 
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Fig. 6. — Observed color-color diagrams at five redshift bins. The total net count is normalized to lxlO 4 in each panel. The intensity 
plot is in counts per pixel of 0.0125x0.0125 mag. The WiggleZ galaxies are plotted as the non-filled contours with a level of 10 galaxies, 
but shifted by A(r f -z f )=l for clarity. The WiggleZ LRR criteria are shown as the black dotted lines, with the r' — z' shifted as well. The 
crosses are the colors of three model galaxies with evolving spectra shown in Fig. [2] and their colors shifted by A (r'-z') =1 are plotted 
as the triangles. Both the WiggleZ galaxies and the neighbors form a continuous sequence at each redshift, with the majority of them 
populated in a blue region. The contamination of dusty star-forming galaxies, indicated by the dashed box, is insignificant for both the 
markers and the neighbors. 



estimate about 5.2±1.0%, 2.4±0.8%, and 0.2±0.5% of 
the galaxies in the three lower-redshift bins in increas- 
ing redshift order may be dusty star- forming galaxies. 
From these fractions, we conclude that dust reddened 
star-forming galaxies are not likely a significant compo- 
nent in the WiggleZ neighbors, at least for z ^0.6. Note 
that the decrease in the dusty galaxy fraction is likely 
due to the shifting of the r' towards the 4000A break 
and possibly the different luminosity depths in the red- 
shift bins, and does not necessarily reflect a real change. 

4.3.2. Control for the Star- Formation Rate of the Markers 

Since the star-formation properties of the WiggleZ 
galaxies themselves are not uniform across redshift due 
to the complex selection criteria, we want to examine 
whether the properties of the neighbors vary with the 
star- formation activity of the WiggleZ markers. We use 
the [OII]A3727A emission line equivalent width (EW) of 
the WiggleZ markers as the proxy for the specific star- 
formation rate of the markers and probe where neighbors 
around markers with different EWQOII] A3727A) popu- 



late the observed color-color diagrams. 

First, we measure the [OIIJA3727A equivalent width of 
the markers. The observed spectra have weak continuum 
due to the short exposure time (~lhr), hence the mea- 
sured EW is noisy for the fainter galaxies. To measure 
the EW([OII]A3727A), we define a window of 10 A cen- 
tered at 3727. 8A as the region for the [OII]A3727Aline. 
The continuum level is determined using a window of 
30 A on each side of the [Oil] line (on the rest frame), 
starting at 3677. 8A and 3757. 8A. Either a linear or 
second-order polynomial function, whichever returns the 
smallest x 2 ? is chosen to describe the fitted continuum 
within the windows. The spectrum is then subtracted 
by this fitted continuum. A bi-Gaussian function is then 
applied based on the data points within all three win- 
dows. The flux of the [OIIJA3727A line is accordingly 
the total net flux under this fitted bi-Gaussian curve. We 
divide the WiggleZ galaxies into three bins based on the 
33.3% percentiles of the EW([OII]A3727A) at each red- 
shift bin. Galaxies without any [OIIJA3727A detection 
are excluded. The median and the la uncertainty of the 
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computed EW([OII]A3727A) are about 112.5±32.6A and 
17.7±8.0A for the highest and lowest EWQOII] A3727A) 
bins at z=0.25-0.35, respectively. 

We re-stack the observed CCM cubes of the neighbors 
by dividing the sample into bins of EWQOII] A3727A) 
within each redshift bin and present them in Figure [7J 
where the WiggleZ markers themselves are again off- 
set by A(r' — z') = 1 for clarity. From Figure [7J it is 
clear that the color-color distributions for neighbors of 
WiggleZ galaxies of different EWQOII]) are very simi- 
lar within the same redshift bin. Kolmogorov-Smirnov 
tests of pairs of both g' — r' and r' — z' distribu- 
tions find no significant difference between the different 
EWQOII] A3727A) bins within the same redshift bin. The 
smallest significant level in all the pair- wise comparisons 
is 0.22. This implies that the properties of the neigh- 
bors do not strongly depend on the properties of the 
WiggleZ galaxies. Furthermore, we find that the N net 
distributions of the neighbors within 0.25 Mpc, although 
not shown here, are identical among all WiggleZ galaxies 
with different EWQOII] A3727A) at a fixed redshift bin. 
This adds to the growing evidence that environment has 
little in fluence in the properties of star-forming galax- 
ies (e.g.. | Balogh et al. 2004; Y ee et al 



2005; Ca rter et al. 
2QQ7HBalogh et al.l 



2001 



iRines et a l. 2005; Cassata et al. 



2009). Therefore, we conclude that the insignificant de- 
pendence between the properties of the neighbors and the 
markers allows us to explore galaxy evolution using the 
neighbors, even though the WiggleZ galaxies may cover 
different ranges of properties at different redshifts due to 
the survey selection criteria; for example, the r / =20-22.5 
criterion naturally select more massive galaxies at higher 
redshifts. 

Since one of the primary WiggleZ target selection crite- 
ria is based on UV flux, we also check whether the NUV 
luminosity of the markers affects the color properties of 
the neighbors. This is essentially testing whether the 
total star-formation rate of the markers affects the color- 
color distributions of the neighboring galaxies. To do so, 
we divide the markers into three groups in each redshift 
bin based on their rest-frame NUV luminosity. Because 
the sample of the markers is not complete to the same 
rest-frame NUV depth, direct comparisons between dif- 
ferent redshift bins cannot be done. However, within the 
same redshift bin, we can compare the color-color distri- 
butions of the markers and their neighbors over its range 
of rest-frame NUV luminosity. The result is presented in 
Fig. El As in the case for EWQOII] A3727A), we observe, 
and confirm with Kolmogorov-Smirnov tests, that within 
a fixed redshift bin the color distributions of the neigh- 
bors around markers with different NUV luminosities 
are statistically identical. This supports our conclusion 
of Fig. [71 that the properties of the neighbors are not 
significantly affected by, or strongly correlated with, the 
characteristics of the markers. We also find that the op- 
tical color distributions of the markers themselves do not 
strongly depend on the absolute NUV luminosity. 

4.3.3. Control for AGN Candidates in the Markers 

Since AGN activities in galaxies are sometime re- 
sponsible for emission lines, we are also interested in 
whether the neighbors of the AGN hosts have similar 
properties as those of normal star- forming galaxies, i.e., 



the rest of the WiggleZ galaxies. The most common 
way to distinguish AGNs and star-forming galaxies is 
based on the ratios of \QIII\/ Hb and [Nlll/iJn, em ission 
lines, the so-called BPT plot ([Baldwin et al.lll981f ). This 
method works only at z < 0.48 for the WiggleZ survey 
where all these emission lines can be de tected within the 
spectr al wavelength coverage. Recently Bongi orno et al.l 
(2010) use a method similar to the BPT plot to separate 
AGNs and star- forming galaxies at 0.50 < z < 0.92, 
i.e., [OIII]A5007A/i^ versus [OH]/ Hp, in the zCOS- 
MOS survey. The separation in this diagnostic dia- 
gram was derived empirically using the observed data by 
studying the positions in the diagram of AGN and star- 
forming g alaxies which wer e classified based on the BPT 
plot (e.g., iRola et!iLlll99ft ILamareille et~alll2004f ). We 
present the [Oil]/ 'Hp vs [Olll]/Hp plot for all our Wig- 



gleZ galaxies at 0.25 < z < 0.75 in Fig. [9l where we over- 
plot the analytical expression of Eq. 3 of [Bon giorno et al.l 
(|2010l ) for the demarcation curves between star-forming 



galaxies and AGNs. We divide the WiggleZ markers into 
three groups based on their locations on the [Oil]/ Hp 
plot. Group A is those located below the analytical ex- 
pression in the star-forming region. Group B contains a 
mix of star- forming and AGN galaxies, located in a nar- 
row strip region centered at the analytical expression. 
Group C is the AGN candidates lying above the ana- 
lytical expression. These groups contribute about 67%, 
20%, and 13% of the galaxies, respectively. 

We find that the neighbors of each group have similar 
N net distribution (Fig. [TO]) , and occupy essentially iden- 
tical regions on the color-color diagram. The N net here 
is comp uted using a limit of M*, + 1 , with M*, derived 
in §4.51 This echoes our conclusion that the properties of 
the neighbors are not strongly affected by the properties 
of the markers themselves. Our results suggest that Wig- 
gleZ galaxies hosting AGN are in environments similar 
to other WiggleZ galaxies. This conclusion is consistent 
with other more detailed studies of Seyfert galaxies that 
the average environment of their hosts is not significantly 
different from other galaxies of simila r properties (e.g., 
iDe Robertis et al.l 11993 iSchmittl 1200 lh . 

4.4. Rest-Frame Color-Magnitude Diagrams 

Figure [11] presents the rest-frame color-magnitude di- 
agrams for the WiggleZ neighbors in redshift bins of 
z ~0.3 to z ~0.7. These are made by summing the el- 
ements of the color-color-magnitude cubes along the x- 
axis (i.e., r'-z' axis). Based on the conclusion in the pre- 
vious two subsections ( §4.3.21 and §4. 3. 3ft that the prop- 
erties of the neighbors are not strongly dependent on the 
properties of the WiggleZ galaxies themselves, the com- 
parisons of the color-magnitude diagrams at difference 
redshifts can provide us useful insights into galaxy evo- 
lution. 

The first observation we can glean from the Figure is 
that the WiggleZ neighbors populate regions of both the 
red sequence and the blue cloud at each redshift bin. 
The majority of them are in the blue cloud. A gap be- 
tween the red sequence and blue cloud is seen. In general, 
the red sequence at each redshift can be approximated 
by a horizontal line of the color of early- type galaxies. 
The flatness of the red-sequence is more likely a result 
of the relatively low signal of the red-sequence galax- 
ies in our data, making them insufficient for deriving an 
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Fig. 7. — Observed color-color diagrams for the WiggleZ and their neighbors where the markers are controlled for their EW([OII] A3727A). 
The total intensity for the neighbors in each panel is normalized to a total count of lxlO 4 , and on the same scale as Fig. [6] The WiggleZ 
galaxies are plotted as the open pink contours with a level of 10 galaxies. The green contours are those in Fig. [6] for all the WiggleZ 
markers at a fixed redshift bin as reference. Within the same redshift bin, both the g' — r' and r' — i' color distributions of the WiggleZ 
galaxies and their neighbors remain similar for samples with different EW([OII] A3727A), indicating that the specific star-formation rate of 
the parent markers does not affect the properties of the neighbors. 

accurate fit, rather than due to the nature of the red 
sequence. The dispersion in the red sequence appears 
to be larger at higher redshift; however, this can be 
mostly attributed to the larger photometric uncertain- 
ties for galaxies in these subsamples. We also plot on 
Figure [TT] the color-magnitude distribution of the Wig- 
gleZ galaxies in each redshift bin as contours, shifted 
by -0.5 mag in g' — r' color. We note that the Wig- 
gleZ galaxies are primarily distributed along the bright 
blue edge of blue-cloud galaxies, reflecting that they are 
strong star- forming galaxies. 



4.5. Galaxy Luminosity Function 

The galaxy luminosity function (GLF) offers a conve- 
nient tool in exploring the different components of the 
galaxy population in a sample. The most wi dely used 
form for the GLF is the Schechter function (S chechterl 
1976), which can be characterized by three parameters: 
the normalization density (/>*, the characteristic magni- 
tude M*, and the faint-end slope a. It has been found 
that a depends strongly on the galaxy SED type. The 
redshift evolution of the GLF, however, is also strongly 



dependent on SED types (e.g . , iWolf et al.ll2Q0l iLiu et al.l 
12008b iSalimbeni eFall 120081 ). The GLF for early-type 
galaxies are described better with a shallower (sometime, 
down- turning) a, and they are more abundant towards 
low redshift. In contrast, late- type galaxies have a GLF 
with a steep a, and their number density is largely un- 
changed toward low redshift. In this subsection, we ex- 
plore the GLF for the WiggleZ neighbors, and investi- 
gate the galaxy population components by examining the 
shape of their GLF and their evolution. 

4.5.1. Constructing the GLF 

The GLF of the neighbors is constructed by project- 
ing the CCM cube along the z-axis (i.e., the r' magni- 
tude) to produce counts as a function of r f magnitude. 
We have conducted and cross-checked the analyses us- 
ing both the observed and rest-frame CCM cubes. Here, 
we present only the results using the rest-frame CCM 
cubes, for which </, r', and z' have been k-corrected to 
z=0.50. We note that the observed r'-band at z = 0.5 is 
approximately equivalent to the rest B band. 

We also separate the neighbors into red and blue pop- 
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Fig. 8. — Similar to Fig. [Tlbut controlling for the markers' NUV luminosity. At a fixed redshift bin, both the colors of the neighbors and 
the markers themselves remain similar regardless of the NUV luminosity. 




log(EW[Oll]X3727A)/EW(Ho) 

Fig. 9. — Selection of AGN candidates. The WiggleZ galaxies 
ar e classified into three g roups by the empirical separation used 
in Bongiorno et al. (2010). The galaxies in group A are primarily 
star- forming galaxies, while those in group C are possibly AGN 
hosts. Group B contains all the rest of the galaxies between group 
A and C, within the regions enclosed by two dashed red lines in 
the plot . 

illations. The division between the red and blue popula- 
tions is chosen to be the g' — r' color half way between the 







Group A 
. Group B 
. Group C 



-20-10 10 20 30 40 -20-10 10 20 30 40 -20-10 10 20 30 40 
N n „/Mpc 2 

Fig. 10. — The N ne t/Mpc 2 distributions of each group in Fig. 
[9] presented as filled blue, thin black, and thick red histograms, 
respectively. The N ne t is computed using the M*, + 1 limit. Each 
histogram has a binsize of 5. All these these groups have similar 
N ne t/Mpc 2 distributions at each redshift. 

non-evolving early-type and the star-forming r model of 
Figure [2] at each redshift, equivalent to 0.27 mag bluer 
than the red-sequence color. To adjust for minor sys- 
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Fig. 11. — The rest-frame color-magnitude of the WiggleZ neighbors at each redshift. The intensity of each panel is normalized to a total 
count of lxlO 4 , and is in units of counts per pixel with a size of 0.0125x0.025 (color x magnitude) mag. The vertical pink dashed lines 
mark the derived °- 5 M*, for the 'All' subsample in ^4.51 The horizontal solid line in each panel marks the separation between blue and red 

galaxies. The color of the red sequence is indicated by the horizontal dotted line. The black vertical dashed lines mark the 5 M r / limit. 
The red sequence of the neighbor galaxies is observable at each redshift. For reference, the WiggleZ galaxies themselves are overplotted as 
the green contours shifted by -0.5 magnitudes in M(g f — r')o.5- 



tematic effects in the photometry, the g' — r' color (at 
rest z = 0.5) of the red-sequence in each redshift bin 
in Fig. [IT] is determined empirically as the reference. 
This is done by examining the g' — r f color distribution 
of galaxies in each redshift bin which are brighter than 
5 M r * + 1 and have colors redder than (g' - r') rs - 0.27, 
where (g f — r') rs is the non-evolving rest-frame model 
early- type galaxy color of Figure [2j The distribution is 
fitted with a Gaussian, and the peak is used as the red- 
sequence color. The red-sequence colors and the bound- 
aries between the red and blue populations are indicated 
in Fig. [Til by the dotted and solid horizontal lines, respec- 
tively. We note that the computed g' — r f red-sequence 
colors are essentially identical to those of the models, 
with the exception of the z = 0.3 bin, where the g' — r' 
separation for red and blue galaxies appears to be ~0.06 
bluer than the computed color of (g f — r') rs — 0.27. 

The GLF results are presented in Figure [12] We denote 
the three subsamples of galaxies and their GLF as 'All', 
'Blue', and 'Red'. The errors in the y-axis in each r' bin is 
computed using Poisson statistics. We also compute the 
GLF using the same method for the RCS-WSF sample 



for comparison with the WiggleZ neighbor galaxies at the 
two lower redshift bins, allowing us to examine the effects 
of environment on galaxy evolution. The CMD and the 
resultant GLFs for the cluster sample are plotted in Fig. 

4.5.2. The Schechter Function Fit 

We fit the GLF discussed above using the Schechter 
function. We determine the completeness limit in 
0-5 M r / for each redshift bin by examining the total net 
counts in 0.1 magnitude bins as a function of 0,5 M r /, 
smoothed by a three-bin kernel. We use the bin 0.1 
mag brighter than the peak as the limit for fitting the 
Schechter function, giving °- 5 M Hm = [-17.9, -18.6, -19.0, 
-19.8, -20.3] for the five redshift bins. These limits are 
marked as vertical purple (dashed) lines in Fig. [12] 

To investigate the shape of the GLF, we first allow a 
to vary in fitting the Schechter function. The results of 
the Schechter function fits are tabulated in Table [3] As 
evident from the CMD (Fig. El, the GLFs in Fig. [12] of 
the WiggleZ neighbor galaxy populations are dominated 
by blue galaxies. It is not surprising that the GLFs of 
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Fig. 12. — The riband luminosity function for all the neighbors of the WiggleZ galaxies. The red, green, and blue curves are the fitted 
Schechter LF with a=-1.50, 1.30, and -1.0, respectively. The Schechter LF for the red neighbors are described by a= -0.70 and -0.40 shown 
as the dotted purple and dashed green curves. The vertical purple line marks the sample complete limit. The GLF for the 'All' and 'Blue' 
neighbors are better described by a steeper a than that to the 'Red' subsample. 



the 'All' and the 'Blue' subsamples are similar, because 
~ 80% of the neighbors are blue galaxies. 

For the 'Blue' and 'All' samples, we find that the fitted 
a appear to become less steep at higher redshift. How- 
ever, this could in part be the result of the different sam- 
pling magnitudes in the subsamples; the shallower abso- 
lute magnitude limits in the higher redshift bins make a 
less well-determined. This is especially true at the high- 
est redshift bin, where the completeness limit is only 
about 1 to 1.5 mag past M*. Discarding the z = 0.7 
bin, the derived a for the 'All' and 'Blue' subsamples 
range between ~ — 1.0 and ~ —1.5, comparable to those 
in th e literature for s t ar- forming; or late- type galaxies 
(e^ iWolf e"taH l200l iChristlein et~aH 120091 : iLiu et al.l 
2008). The possibility that the apparent changes in a 
at different redshifts are due to the change in the sam- 
pling limit can be demonstrated by fitting the three low- 
redshift bins for the 'Blue' subsample to a common limit 
of 5 M r / = —19.0, which is complete for the three bins. 
We obtain the results: a = [-1.15 ± 0.17, -1.38 ± 0.16, 
-1.03 ±0.12] for z = [0.3, 0.4, 0.5], which are consistent 
with being identical at ~ 1.75<r. 

We overplot the Schechter LF in Fig. [12] with fixed 
a=[-1.50, -1.30, -1.00] fitted for the 'Blue' subsamples 



to show the ability of the data to distinguish between 
the faint-end slopes within this range. By combining the 
three lower-redshift bins (0.25 < z < 0.55) where the 
data are of sufficient depth to obtain a robust measure- 
ment of a, we find for the 'Blue' galaxies a best-fit a of 
— 1.18 ± 0.08. For consistency of comparison of 0,5 M r / 
among redshift bins, we adopt an a of —1.3 for all red- 
shift bins, and tabulate the fitted °- 5 M*, in Table 03 

The best fitting results for °- 5 M*, and a for the Wig- 
gleZ 'Red' neighbor galaxy samples at different redshift 
bins are also presented in Table 03 For the lower redshift 
bins, where the data are of sufficient depth, the GLF is 
considerably less steep at the faint-end compared to that 
of the 'Blue' subsample. At higher redshifts, the red 
GLFs can be described by a steeper a, but this is again 
likely because the faint-end of the GLF is not observable 
at these redshifts. We use the three lower redshift bins 
to establish a more robust estimate of the faint-end slope 
of the red galaxy GLF. We obtain a = -0.45 ± 0.13 by 
summing the data (to 5 M r / = —19.0) in these bins. 

For the purpose of comparison, we perform the same 
analysis for the RCS-WSF sample of markers in RCS2 
clusters. Here, the CMD shows a very strong red se- 
quence and the total GLF is dominated by red galaxies, 
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Fig. 13. — The color-magnitude diagrams (Top) and luminosity 
functions (Bottom) for the RCS-WSF sample for the lowest two 
redshift bins. The red-sequence galaxies are domina nt in both the 
plots. The symbols are of the same meaning as Fig. [12] 

as shown in Figure [T3]) . The faint-ends of the 'Blue' 
galaxy GLF in the two redshift bins have a similar a as 
those for the corresponding subsamples of the WiggleZ 
neighbors. For the purpose of comparing M*, we also 
fit the RCS-WSF 'Blue' samples using a = [-1.50, -1.30, 
-1.00] with the results listed in Table H 

Opposite to the WiggleZ neighbor samples, the GLF of 
the RCS-WSF sample of cluster galaxy targets is dom- 
inated by the 'Red' galaxy sample. The faint-end slope 
for these red galaxies appears to be significantly steeper 
than that of the WiggleZ counterpart, with a steeper 
than -0.65. The combined 'Red' GLF for the two red- 
shift bins of the RCS-WSF data produces a best fitting 
a = —0.74 ±0.06. A similar combination for the WiggleZ 
neighbor sample produces a — —0.18 ± 0.17. In order to 
provide a consistent comparison for the 0,5 M*, values for 
the 'Red' WiggleZ neighbor sample and the 'Red' RCS2 
cluster neighbor sample, we also refit all the red galaxy 
subsamples using a = —0.4 and -0.7, and the results are 
tabulated in Table [3] 

4.5.3. The Evolution of 5 M*, 

We examine the evolution of the M*, parameter of the 
Schechter LF by assuming a simple li near dependence 
between M*, and redshift, as used by iLin et all (jl999) 
and others. We can write: °- 5 M*, (z)=m*, — DM — (k — 
fc^o.50) + Q(z — 0.5), where DM is the distance modulus 
and k is the k-correction. The evolution term Q is de- 
rived from a linear fit to 0,5 M* f as a function of redshift 
in the form of °- 5 M* f (z) = °- 6 M*,(z = 0.5) + Q(z - 0.5). 
The derived Q depends on the value of a; a steeper 
a results in a smaller Q. Because the fitted a and 
M*, correlate with each other and a are marginally 
different in our redshift divisions, we use the M*, de- 
rived with a = —1.30 for the 'All' and 'Blue' subsam- 
ples and a= -0.40 for the 'Red' neighbors, obtaining 
Q=[-1.31±0.24 ,-2.10±0.26 ,-1.59±0.30] for the 'All', 
'Blue', and 'Red' subsamples, respectively. The results 
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Fig. 14. — ' 5 M*, as a function of redshift derived from Fig. [12] 
using three fixed ol. The dashed curve with triangles is the Mg 
from Faber et al. (2007), and the one with crosses symbols is the 
5 M r / by applying a simple color transfer from B to r' . The linear 
fits for the WiggleZ neighbor galaxies are overplotted as the dotted 
lines. The RCS-WSF sample is overplotted as purple open circle. 
The derived evolutionary term, i.e., the slope of the fit, is -1.31 for 
all the neighbors (the left panel), -2.10 with ol— -1.30 for the blue 
neighbors (the middle panel), and -1.59 for the red neighbors with 
a=-0.40. 

are shown in Fig. [14] and tabulated in Table [3] We also 
derive Q using a = — 1.0 and -1.5 for the 'All' and 'Blue' 
neighbors, and a = —0.7 and -0.4 for the 'Red' subsam- 
ple. These fits are shown in Fig. [14] as dotted line for 
each redshift bins. 

We adopt the parametrization of °- 6 M*,(z) = -21.44- 
1.31 — 0.5) from the 'All' sample in the derivation of 
the red-galaxy fraction in the following subsection. 

4.6. The Fraction of Red Neighbors 

In Figure [6] we have observed that dusty star- forming 
galaxies are not common for both the WiggleZ galaxies 
and their neighbors. The red neighbors of the WiggleZ 
galaxies are more likely to be those which have completed 
their star formation. We investigate the fraction of these 
red passive neighbors (f re d) as a function of redshift. We 
define red neighbor galaxies the same way as described in 
£ 34.5.11 The boundaries separating the blue and the red 
galaxies are shown in Fig. [TTJ which is about 0.27 mag 
bluer than the red-sequence, k-corrected to z = 0.5. This 
color separation is equivalent to the gap in the galaxy 
bimodal color distribution at a fixed redshift seen in our 
samples. The fraction of net red neighbors to the total 
net neighbor counts as a function of redshift is presented 
in Fig. [15] We compute f re d to four different depths 
of 5 M* + AM with 5 M r * as described in 34X31 and 
AM=[0.0, 0.5, 1.0, 1.4], with the largest AM determined 
by the depth of sampling of the data from the largest 
redshift bin. The error in f rec i is estimated using Poisson 
statistics. 

We observe that f re d is not a strong function of red- 
shift. Within the uncertainties of the measurements, 
f re d can be described as flat between z = 0.25 and 0.65, 
with an average of ~0.28±0.01 for the magnitude limit 
of °- 5 M r *, + 1.4 . At z ~ 0.7, there appears to be a drop 
in f re d at ~ 3.5<j level. We also see a relatively small 
change in f re d within the relatively magnitude range we 
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Fig. 15. — Left: Fraction of red neighbors f re d as a function of redshift to different M r i de pth. Differen t curv es represent the f re d 
computed using different magnitude depth as noted in the plot. The gray straight lines are from Iovino et al. (2010) for their Sample III 
(< M* + 0.8) 'All' (solid linestyle) and 'Isolated' (dashed linestyle) subsamples. The f red of the RCS-WSF sample at z ~0.30 and 0.40 
are overplotted as individual open symbols, which have a value of ~0.8. The f re d for the z = 0.75 — 0.85, which is beyond the redshift 
range of our main sample, are plotted as individual points for magnitude limits of M*, and M*, +1.2 (black square) to illustrate the 

continuous decrease of f re d beyond z ~ 0.7. Right: f re d computed using °- 5 M*, + 1 depth after controlling for the EW([OII] A3727A) of 
the WiggleZ markers. The gray shaded region is the f re d using all markers with the same magnitude limit. The f re d is similar regardless 
of the properties of the markers. 



probed - changing by the order 0.1 over the 1.4 mag- 
nitude range, often within the uncertainties of the mea- 
surements in the same redshift bin. However, we note 
that differentially, as indicated by the GLFs of the red 
and blue galaxies, the fraction of blue galaxies increases 
significantly at the faint-end. Using data combining the 
two lowest redshift bins where we can sample down to 
°- 5 M* + 2.5, we find f red ~ 0.22 ± 0.01, a difference of 
- 5a from the f red ~ 0.45 ± 0.04 measured at 5 M r * . 

For comparison, the f re d for galaxies around the RCS- 
WSF sample of cluster galaxy markers are also plotted 
in Fig. [15j which have an average f re d ~ 0.8 for the 
z ~0.3 and 0.4 bins. As expected, the f re d for the RCS- 
WSF sample is much larger than that for the WiggleZ 
neighbors. 

We also compute the f re d values for subsamples of Wig- 
gleZ markers in different EW([OII]A3727A) bins. We di- 
vide the markers into three groups based on the 33.3 and 
66.7 percentiles in the distribution of EWQOII] A3727A). 
The f re d values, computed using the limit of °' 5 M*, + 1, 



as a function of redshift are presented in the right 
panel of Figure. [151 The gray shaded area is the 
f re d calculated with the same magnitude limit using all 
the neighbors (i.e., the dashed blue curve in the left 
panel), overlaid for comparison. Although the Wig- 
gleZ markers with stronger EW([OII]A3727A) appear 
to have somewhat lower f re d, the f re d for the different 
bins are well within the individual measurement uncer- 
tainties, with the exception of the z = 0.5 bin. Av- 
eraging over all redshift bins, we find the mean f re d 
for the three EW([OII]A3727A) bins, from low to high, 
to be 0.305 ± 0.019, 0.278 ± 0.021, and 0.246 ± 0.021, 
or, a difference of about 2a between the weakest and 
strongest [OIIJA3727A samples. Thus, there is some ev- 
idence that the neighbors of markers with the largest 
EW([OII]A3727A may have a lower average f re d, indi- 
cating regions around galaxies with large specific star- 
formation rates may have a larger fraction of star- 
formation galaxies. However, the relatively low signifi- 
cance difference (which comes mostly from the z = 0.5 
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TABLE 3 

Parameters for Luminosity Functions of the WiggleZ 
Neighbors 



All Neighbors 


redshift 




a 




fitting a 








0.25-0.35 


-21.46±0.18 


-1.42±0.05 




0.35-0.45 


-21.33±0.21 


-1.35±0.08 




0.45-0.55 


-20.89±0.14 


-0.89±0.10 




0.55-0.65 


-21.23±0.15 


-0.97±0.12 




0.65-0.75 


-21.10±0.16 


-0.76±0.18 




redshift 


°- 5 M*, 


0.5 M * 


0.5 M * 
r' 




a= -1.50 


a= -1.30 


a= -1.00 


0.25-0.35 


-21.68±0.11 


-21.16±0.08 


-20.60±0.06 


0.35-0.45 


-21.70±0.12 


-21.22±0.09 


-20.69±0.06 


0.45-0.55 


-21.92±0.10 


-21.51±0.08 


-21.02±0.05 


0.55-0.65 


-21.98±0.09 


-21.65±0.07 


-21.26±0.05 


0.65-0.75 


-21.89±0.08 


-21.63±0.06 


-21.30±0.05 


Q 


-0.64±0.32 


-1.31±0.24 


-1.95±0.19 


°- 5 M*,(z = 0.5) 


-21.84±0.05 


-21.44±0.04 


-20.98±0.03 


Blue Neighbors 


redshift 


°- 5 M*, 

r' 


a 




fitting a 








0.25-0.35 


-21.41±0.25 


-1.56±0.06 




0.35-0.45 


-21.34±0.29 


-1.49±0.10 




0.45-0.55 


-20.90±0.17 


-1.02±0.11 




0.55-0.65 


-21.19±0.19 


-1.05±0.15 




0.65-0.75 


-21.02±0.18 


-0.74±0.21 




redshift 


o.5 M * 


0.5 M * 


°- 5 M*, 




a= -1.50 


a= -1.30 


a= -1.00 


0.25-0.35 


-21.23±0.11 


-20.76±0.08 


-20.23±0.06 


0.35-0.45 


-21.34±0.13 


-20.90±0.10 


-20.40±0.07 


0.45-0.55 


-21.72±0.11 


-21.32±0.09 


-20.87±0.06 


0.55-0.65 


-21.80±0.10 


-21.49±0.08 


-21.13±0.06 


0.65-0.75 


-21.79±0.09 


-21.55±0.07 


-21.24±0.06 


Q 


-1.49±0.34 


-2.10±0.26 


-2.71±0.20 


°- 5 m*,(z = o.5) 


-21.58±0.05 


-21.21±0.04 


-20.78±0.03 


Red Neighbors 


redshift 


0.5 M * 
r' 


a 




fitting a 


0.25-0.35 


-20.23±0.20 


0.099±0.27 




0.35-0.45 


-20.74±0.22 


-0.50±0.19 




0.45-0.55 


-20.60±0.20 


-0.26±0.21 




0.55-0.65 


-21.21±0.24 


-0.64±0.22 




0.65-0.75 


-21.33±0.38 


-0.77±0.37 




redshift 


°- 5 M*, 
a= -0.70 


a= -0.40 




0.25-0.35 


-20.51±0.12 


-20.40±0.09 




0.35-0.45 


-20.96±0.10 


-20.63±0.09 




0.45-0.55 


-20.99±0.09 


-20.73±0.08 




0.55-0.65 


-21.28±0.09 


-20.98±0.07 




0.65-0.75 


-21.26±0.11 


-21.01±0.09 




Q 


-1.79±0.37 


-1.59±0.30 




°- 5 M*,(z = 0.5) 


-21.02±0.05 


-20.76±0.04 





bin) is consistent with our conclusion of Figures and 
[8] that the properties of the neighbors are not strongly 
dependent on the characteristics of the WiggleZ galaxies 
themselves. 



5. DISCUSSION 

5.1. Robustness of Our Results 

We have demonstrated that our novel approach - prob- 
ing photometric properties of galaxies using color-color- 
magnitude cubes of neighbors around galaxies of known 



TABLE 4 

Parameters for Luminosity Functions of the RCS-WSF 
Sample 



All Neighbors 


redshift 


°- 5 M*, 
r' 


a 




fitting a 








0.25-0.35 


-21.04±0.08 


-1.10±0.03 




0.35-0.45 


-21.05±0.10 


-0.94±0.06 




redshift 


0.5 M * 


o.5 M * 

r' 


o.5 M * 

r' 




a= -1.50 


a= -1.30 


a= -1.00 


0.25-0.35 


-22.14±0.07 


-21.53±0.05 


-20.81±0.03 


0.35-0.45 


-22.14±0.09 


-21.67±0.06 


-21.13±0.05 


Blue Neighbors 


redshift 


°- 5 M*, 

r' 


a 




fitting a 








0.25-0.35 


-20.94±0.41 


-1.52±0.03 




0.35-0.45 


-21.59±0.73 


-1.55±0.06 




redshift 






- 5 M*, 




a= -1.50 


a= -1.30 


a= -1.00 


0.25-0.35 


-20.87±0.18 


-20.30±0.13 


-20.00±-0.0 


0.35-0.45 


-21.39±0.29 


-20.82±0.21 


-20.20±0.15 


Red Neighbors 








redshift 


o.5 M * 
r 


a 




fitting a 








0.25-0.35 


-20.52±0.09 


-0.64±0.06 




0.35-0.45 


-20.81±0.09 


-0.66±0.06 




redshift 


0.5 M * 
a= -0.70 


a= -0.40 




0.25-0.35 


-20.62±0.04 


-20.14±0.04 




0.35-0.45 


-20.86±0.04 


-20.45±0.03 





redshift - has yielded robust and interesting results. The 
advantages of our method are: (1) being straightforward, 
as it is equivalent to counting galaxies within an aperture 
and applying statistical background corrections; (2) not 
sensitive to whether the sample of the markers is com- 
plete; (3) allowing us to derive the photometric proper- 
ties of galaxies to a limit considerably deeper than the 
corresponding limits for the spectroscopic sample; and 
(4) providing a complete census of the neighboring galax- 
ies independent of the spectral/color properties of the 
galaxies. The key point in the method is the assumption 
that the marker galaxies and the excess galaxy counts 
around them are at the same redshift. The net neighbor 
counts around an individual galaxy are not statistically 
meaningful, but stacking color-color-magnitude cubes of 
a large number of markers provide statistically significant 
quantities. However, much care is required in carrying 
out the procedure. 

First, proper background correction is essential in our 
method, especially when the signal of the net excess is 
only a fraction of the total galaxy counts. Any system- 
atic discrepancy in the background subtraction would 
have a profound effect on the result. If the background 
is under-subtracted, the intrinsic properties of the net 
excess galaxies will be overwhelmed by the background 
counts; for instance, it may result in a power-law galaxy 
luminosity function without an apparent bend/knee. On 
the other hand, over-subtracting the background may 
result in small, or even negative, net counts, and hence 
no analyses can be done. The very large angular size 
of the RCS2 patches allows us to use photometric data 
from the same patch as the markers to estimate the back- 
ground correction, minimizing any possible systematic is- 
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sues. The use of uniform random catalogs to map out the 
imaging survey area, CCD chip gaps, bad CCD columns, 
bright star halos and other artifacts, ensure the proper 
treatment of the sampling aperture size. Our exercise 
of measuring excess cou nts centered on a large number 
of random positions in §2.21 demonstrates that we have 
handled the background correction properly, as the net 
excess around the randomly drawn points has a mean 
value of essentially zero and is not a function of galaxy 
magnitudes or colors. 

Another way of verifying the robustness of our back- 
ground subtraction technique is to use the WiggleZ red- 
shift sample. While it poorly samples the whole pho- 
tometric galaxy catalog, the very large sample of Wig- 
gleZ redshifts allows us to test the principle of back- 
ground subtraction in general, and our technique in par- 
ticular. We perform this test by comparing the ratio of 
the net excess galaxy counts to the total counts within 
the r p = 0.25Mpc aperture in both the photometric and 
WiggleZ redshift samples. We note that this method may 
not produce an exact comparison, as the WiggleZ sample 
has a number of selection criteria which may produce dif- 
ferent selection effects for galaxies at different redshifts 
that are not possible to mimic using a purely magnitude- 
limited photometric sample. The most significant se- 
lection effect that produces a significant difference in 
the redshift distribution between a r'-band selected sam- 
ple and the WiggleZ sample is the low-redshift-rejection 
(LRR) applied to the selection of WiggleZ targets (see 
§2.3). Thus, we limit our comparison to using data from 
WiggleZ markers in the three high-redshift bins of our 
sample (0.45 < z < 0.75). For the photometric sample 
of counts of neigbors to the WiggleZ markers, we count 
only galaxies with 20 < r' < 22.5 (which is the Wig- 
gleZ optical-band selection criterion). For the WiggleZ 
galaxies, we count all WiggleZ glaxies within the aper- 
ture of each marker, and deem galaxies with a redshift 
within 0.002(1 + z) of a marker galaxy as associated. We 
find the ratio of net excess counts to total counts in the 
aperture to be 0.193 ±0.010 for the photometric sample, 
where the uncertainty is based on Poisson statistics. For 
the equivalent ratio from the WiggleZ redshift sample, 
we obtain 0.196 ± 0.014 (with 235 out of 1202 galaxies 
satisfying the redshift criterion). Thus, our statistical 
background subtraction technique produces net counts 
that are entirely consistent with results using a sample 
of galaxies with known redshifts, further adding confi- 
dence to the robustness of method. 

Second, incorrect redshift measurements of the mark- 
ers will dilute the results. The redshift plays the key role 
in computing the distance modulus, k-correction, and the 
aperture size for a fixed physical diameter. The first two 
impact directly on the rest-frame color and magnitude 
distributions of the neighbors. The latter affects the sur- 
face density of the excess counts. Since the redshift confi- 
dence for the markers with redshift quality flag Q ZS pec=3 
is ~78.7%, we test the effect of incorrect redshift mea- 
surement by repeating our analyses using objects with 
Qzspec =/ ^ and 5 only, at the cost of having a weaker sig- 
nal due to stacking fewer color-color-magnitude cubes. 
We find that the results are very similar to what we have 
presented here. 

Finally, an interesting question is how often we have 
another WiggleZ galaxy within the aperture (r p = 0.25 



Mpc) centered at a WiggleZ marker. We find that ~5- 
10% of the markers have at least one other WiggleZ 
galaxy within a radius of 0.25 Mpc, but the percent- 
age drops to ~1% if these other WiggleZ galaxies have 
to be at a similar redshift (Az < 0.002(1 ± z)) as the 
marker. The range in the percentage is patch depen- 
dent, as the number density of WiggleZ galaxies varies 
in different patches. We also find that the number of 
the enclosed WiggleZ galaxies is not a strong function of 
redshift. Therefore, we conclude that the properties of 
the WiggleZ neighbors should not be significantly con- 
taminated by the WiggleZ galaxies themselves. 

5.2. The WiggleZ Neighbor Galaxy Luminosity Function 
5.2.1. The Schechter Function Fit 

Our WiggleZ galaxy neighbor sample represents a com- 
plete census of galaxies in the neighborhood of star- 
forming galaxies over the redshift range of 0.25 to 0.75. 
On the color- magnitude plane, we separate the galaxy 
sample into red and blue galaxy subsamples. The GLFs 
for both the 'Blue' and 'Red' galaxy samples can be fit- 
ted very well with single Schechter functions. Based on 
the low-redshift bins, where there is sufficient depth to 
measure the faint-end slope unambiguously, they have 
different shapes such that the red galaxies have a dipping 
faint-end (best fitted with — 0.6), whereas the blue 
galaxies have a steep rising faint-end which is best fitted 
with a ~ —1.3. This is similar in general to what is seen 
in clusters, where red-sequence galaxies typically have a 
turn over in the GLF, while the blue ga laxies increase in 
number steeply at the faint end (e.g., iBarkhouse et al.l 
[20071 : IWolf et^l2009l) . 

We can make a direct comparison between the Wig- 
gleZ neighbor sample, representing galaxies associated 
with star- forming galaxies, and the neighbors of the RCS- 
WSF sample, representing galaxies in dense cluster envi- 
ronments, by combining the z ~ 0.3 and 0.4 subsam- 
ples. We find the 'Blue' galaxy samples to have es- 
sentially identical GLF parameters: a = — 1.42 ± 0.08 
and -1.57 ± 0.13, and 5 M r * = -21.15 ± 0.19 and 
-21.37 ±0.47. The red galaxy GLFs for the two samples 
appear to have different faint-end slopes at the 4.3<r level, 
with a = -0.18 ±0.17, and -0.74 ±0.06 for the WiggleZ 
neighbors and the RCS-WSF neighbors. The RCS-WSF 
red galaxy neighbors also have a marginally brighter 
5 M r * of -20.78 ± 0.09, compared to -20.42 ± 0.15 for 
the WiggleZ neighbor sample, at the 2a level. However, 
this difference is most likely a reflection of the different 
a used in the best fits; fitting both 'Red' GLFs using 
a common a = —0.5 for z = 0.25 to 0.45, we obtained 
5 M* =-20.42 ± 0.09 and -20.53 ± 0.07 for the RCS- 
WSF and the WiggleZ neighbor sample, respectively. We 
will further discuss the difference in the GLF shape in 

5.2.2. The Evolution of the Galaxy Luminosity Function 

We use the parameter Q (see ^4.5j) to measure the evo- 
lution of M*, with the asssumption of a being constant 
with redshift for the faint-end of the GLF. The results 
for different samples and a's are shown in Fig. [TH There 
is a general brightening of M* with larger redshift. An 
intereting trend is that that red-sequence galaxies may 
have a lower Q value than that of blue-cloud galaxies, 
~ — 1.6 versus ~ —2.1, but only at the 1.3a. 
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Using a sample drawn from the DEEP2 and COMBO- 
17 data and with SDSS data as the local universe epoch, 
iFaber et all (|2007f ) studied the evolution of the GLF in 
the redshift range of 0.1 < z < 1 in the Mb band, which 
is similar in rest wavelength to our 0,5 M r / band. They ob- 
tained Q values of -1. 23±0.36, -1.34±0.22, -1.20±0.21, 
over the range of z ~ to 1, for their 'All', 'Blue', and 
'Red' samples. These values appear to be lower from 
those derived in our samples (see Table [3]) at moderately 
significant levels, especially for the 'Blue' sample. How- 
ever, the comparison is much more similar if we limit 
their data to the same redshift range as ours. In Fig- 
ure [lU we plot their M B data points. For a more direct 
comparison, we also convert their B band to our 0,5 M r / 
by applyi ng color corrections base d on B — r' colors from 
GISSEL (jBruzual fc Charlotll2QQ3l ) to convert M B to M r > 
at zero redshift, and then k-correct M r > to 0,5 M r /. These 
data are also plotted in Figure [14} 

We refit the Q factor from lFaber et al.l (12007) us- 
ing their data within the redshift range covered by our 
the WiggleZ neighbor sample. We find Q = [—1.45 ± 
0.25, -2.02 ± 0.32, -0.67 ± 0.38] (in their B band) for 
their 'All', 'Blue', and 'Red' samples, respectively, com- 
pared to our values of -1.31 ± 0.24, -2.10 ± 0.26 (a=- 
1.3) and -1.59 ± 0.30 (a = -0.4) in Table O Thus, the 
derived Q values from the two studies over the same red- 
shift range are similar, especially for the 'All' and 'Blue' 
samples. For the 'Red' samples, the WiggleZ neighbors 
have a steeper evolution, at the ~ 2a level. We note that 
over this redshift range the Faber et al. results also pro- 
duce a relatively significant lower Q value for their 'Red' 
sample compared to that for their 'Blue' sample, at the 
2.7<j level, reinforcing a possible similar trend in the Wig- 
gleZ samples. The marginally steeper evolution of the 
'Red' samples in the two data sets could be a reflection 
of the difference in the selection for the red galaxy sam- 
ples. This possible discrepancy could be an indication 
that there is a difference in the evolution of early-type 
galaxies based on the environment. The red galaxies in 
the WiggleZ neighbor sample are primarily in low den- 
sity regions, conducive to star-formation; while those in 
the Faber et al. sample include red galaxies in all en- 
vironments, with an expected bias towards high density 
regions. Thus, it is perhaps not surprising that the two 
samples have different Q values, with the red GLF in low- 
density regions evolving more rapidly. However, we note 
that the differences discussed in this subsection are at 
the ~ 2a levels. Considerably more detailed studies are 
need to firmly establish the differences in the evolution 
of the GLF of different galaxy populations in different 
environments. 

5.3. The Red- Galaxy Fraction 

When the star forming activity in a galaxy ceases, the 
galaxy ought to become redder in colors. In Fig. [T5l we 
explore the redshift dependence of the red galaxy frac- 
tion, f re d. Since the dependence of f re d on z is similar 
for samples of different depths, for the remainder of the 
analysis, we use the °- 5 M r /< 0,5 M*, + 1.4 sample, which 
has f re d measurements with the smallest error bars and is 
complete for all the redshift bins, unless noted otherwise. 

5.3.1. The Redshift Dependence of f re d 



We observe that f re d is remarkably similar at 0.25 < 
z < 0.65, with a value of ~ 0.28. A linear fit to the 4 
points gives a slope of —0.070 ± 0.098. The f re d drops 
to about 0.20 for the z = 0.7 bin. This drop, between 
the z = 0.6 and 0.7 bins, is statistically significant at the 
~ 3.5<j level. If we assume the best linear fit for f re d from 
the four data points with z < 0.65, the f re d at z = 0.7 
is 4.8<j below the extrapolation from the fit. An alterna- 
tive simpler description of the trend is a linear decrease; 
however, this linear fit, to all 5 points, has a reduced x 2 
of 3.1, indicating that it is not a good description of the 
data. 

To investigate whether there is a change in the depen- 
dence of f re d on z at z ~ 0.7, we extend the measurement 
of fred to z ~ 0.8. The WiggleZ sample has a large num- 
ber of galaxies at 0.75 < z < 0.85, which we have not 
included in our analysis because the RCS2 photometry 
is only complete to a relatively shallow absolute magni- 
tude limit of -20.6. Nevertheless, this bin is complete 
to 0,5 M r /+1.2, and hence, we perform the same analysis 
using the 5819 WiggleZ galaxies in this redshift bin as 
markers. The f re d computed using this sample is plot- 
ted on Figure [T5l It shows a continuous drop from the 
z = 0.7 data point. A linear fit applied to all 6 data 
points from z = 0.3toz = 0.8 with a limiting magnitude 
of M* + 1.2, produces a reduced x 2 of 6.3, indicating a 
poor fit. Thus, the addition of the higher redshift data 
adds confidence to the conclusion there is an on-set of a 
significant drop in the red-gala xy fraction a t z > 0.65. 

We compare our results to Iovino et al.l (j2010h . who 
derived the blue galaxy fraction Fu U e for galaxy samples 
in different environments from the zCOSMOS survey. Of 
particular interest is their 'All' and 'Isolated' subsamples 
of Sample III, with a depth of M = M* + 1.5, c o vering 
the redshift range of - 0.2 to 0.6. Ilovino et al.l ([201(1 ) 
parametrize the evolution of Fu U e by a power law of the 
form equivalent to: f re d(z) = 1 — Fbiue(Q)0- + Z Y \ with 
Fbiue(0) = 0.58 ± 0.10 and /3 = 0.69 ± 0.44 for the Iso- 
lated sample, using 3 data points. We plot their Fbiue(z) 
as fred on Figure [15] using their parametrization. Our 
fred has values more similar to their Isolated sample, and 
well below their 'All' sample. This is consistent with our 
expectation that the WiggleZ neighborhoods have by- 
and-large a low galaxy density environment, dominated 
by star- forming galaxies. 

For a more direct comparison, we fit our f re d with the 
powe r law over a l o nger r edshift range (z = 0.20 — 0.70) 
than Ilovino et al.l (|2010D . and obtain (for M* + 1.4): 
F b iue(0) = 0.59 ± 0.04, and f3 = 0.52 ± 0.13 with a re- 
duced x 2 ~ 3.1. Thus, our data show a similar, rate of 
dec rease of f rP A with increasing redshift as that found 
by Ilovino et al.l (j2010f ). However, we note that both 
of the fitting models to our data-the linear fit and the 
power-law fit-have reduced x 2 s considerably larger than 
1: ~ 3 — 4 for both using data up to z = 0.7 and 
z = 0.8. This indicates that a simple continuous decrease 
of fred with redshift is likely not a good description of 
its evolution. 

Thus, an interesting result is that ou r data, cov- 
ering a longer redshift baseline than the Ilovino et al.l 
(|2010f ) study, show a more or less constant f re d up to 
z ~ 0.6 before seeing a drop. Such a description of 
the change in f re d for galaxies in poor environments 
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with redshift is in fact also consisten t with the data of 
the 'Isolated' sample of llovino et al.l (j201Q[ ). with their 
three F^iue data points covering the redshift between 
0.2 and 0.6 being consistent with having similar val- 
ues within their uncertainties. We note that the more 
or less constant f re d~ 0.3 for z<0.6 for the WiggleZ 
neighbor samples is similar to that obtained for galax- 
ies i n low- density r egions at the local universe at z ~ 0; 
e.g., iBalogh et all ([2004( 1 derived f red ~ 0.35 from the 
SDSS sample for their low local-galaxy density samples 
of galaxies of My ^ — 20. The f re d for the WiggleZ galaxy 
neighbors extrapolate to ~ 0.32 at z = 0. Thus, there 
appears to be a relatively small amount of evolution in 
f re d in low galaxy density regions from z ~ 0.6 all the 
way to z ~ 0. 

Because the WiggleZ sample is selected in part by ob- 
served UV and optical fluxes of the objects, this abrupt 
increase in the change in f re d between the 0.6 and 0.7 
redshift bins could conceivably be contributed by some 
unknown correlation or thresholding effect between the 
star-formation rate or luminosity of the markers and 
the properties of their neighbors. To test the effect of 
the UV flux selection criterion, we select two subsam- 
ples of markers within an identical UV absolute lumi- 
nosity range which is complete in the both z ~ 0.6 and 
z ~ 0.7 bins {M^uv < —20.7), and compute their f re d- 
We obtain / red =0.293±0.022, and 0.207±0.017 for the 
neighbors in the two redshift bins, indicating a similar 
and significant (at ~ 3a) drop to that when the whole 
sample without control on M^uv is used. Similarly, we 
compare the f re d of neighbors of subsamples of markers 
with the same r f absolute magnitude ( 0,5 M r / between 
-21.5 and -22.5) in the last two redshift bins, and ob- 
tain / re ^=0.255±0.045 and 0.176±0.033; again, showing 
a drop similar to that obtained using the whole redshift 
bins. Thus, we can conclude that the drop in f re d be- 
tween z ~ 0.6 and 0.7 is likely not a result of the different 
star-formation properties of the central markers. 

5.3.2. The Evolution of the Red- Galaxy Fraction 

Our data allow us to look at the redshift evolution of 
the red-galaxy fraction f re d up to z ~ 0.7 using samples 
complete to an absolute luminosity about 1.5 mag be- 
yond M*. Under the simplest assumption of a sample of 
galaxies in a closed volume, f re d allows us to follow the 
end result of the process of galaxies having their star for- 
mation quenched and eventually tur ning red, becoming 
a me mber of the red sequence (e.g., iKodama fc Bowerl 
2001). However, the use of a luminosity limit produces 
ambiguities in how to interpret the average change in 
the galaxy population between the different epochs, since 
galaxies of a given mass may enter and leave the sam- 
ple depending on their star- formation state. To assist 
in the interpretation of the observed change in f re d with 
redshift, we illustrate in Figure [16] the possible flows of 
galaxies into and out of luminosity-limited blue and red 
galaxy samples. 

In a simple picture, the number of galaxies in the blue 
cloud can be augmented by new star formation in galax- 
ies fainter than the luminosity limit, boosting them into 
the sample, adding to the blue galaxy counts. This re- 
plenishment is represented by the white arrows (A) in 
Figure [TH In a field situation, where infall in general is 
not a major process, one can consider this replenishment 
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Fig. 16. — Schematic color- magnitude diagram showing the flows 
of galaxies into and out of the blue cloud. The narrow red ellipse 
represents the location of the red sequence; the blue ellipse marks 
the blue cloud. The vertical dotted line indicates the magnitude 
limit of the sample for computing f re d- The white arrows (A) show 
the replenishment of galaxies in the blue cloud due to an increase in 
the star-formation rate of blue galaxies fainter than the magnitude 
limit. The black three arrows illustrate the possible flows of galax- 
ies from the blue cloud into the red-sequence after star formation 
quenched. Arrow B3 represents star-formation being quenched by 
an approximately equal-mass merger. Arrow B2 stands for the 
quenching of star formation, while arrow Bl represents the case 
when the galaxy moves back beyond the magnitude limit again 
due to its quenched star formation. 

generally being controlled by the canonical down-sizing 
scenario of star formation - lower mass galaxies may start 
forming stars, and in some instances, boosting the lumi- 
nosity of the galaxy into the blue cloud sample. 

On the other hand, one would expect a continuous flow 
of galaxies from the blue cloud into the red sequence as 
galaxies evolve. This would primarily be created by the 
quenching of star formation in galaxies in the blue cloud 
by various processes. There are several paths that this 
transformation into the red sequence may take, and the 
are illustrated by the black three arrows (Bl, B2, and B3) 
schematically. The luminosity of a galaxy will generally 
fade as its star formation is quenched and turns red (Bl 
and B2), and some of them will become fainter than the 
sample magnitude limit and leave the sample. There 
is also a possibility that the galaxy becomes brighter if 
the cessation of star formation follows a major merger 
(arrow B3). However, in poor environ ments, it is unlikely 
that this is a dominant event; e.g., Hsie h et~a"T1 ([2008) 
found approximately 6% of galaxies may have undergone 
a major merger since z ~ 0.8. 

In the WiggleZ sample of neighbors of strong star- 
forming galaxies we find that f re d is essentially flat from 
z ~ 0.6 to the present, with evidence of a decrease at 
z ^0.6. In the simplest picture, the change in f re d over 
redshift can be interpreted as a change in the relative flow 
rates of galaxies into the red sequence and into the blue 
cloud. Thus, at z^0.6 we can effectively conclude that 
the rate of transformation of blue galaxies into red galax- 
ies, which builds up the red-sequence, is approximately 
equal to the replenishment rate of blue galaxies brighter 
than the sampling limit; while at z ^0.6 the rate of trans- 
formation of blue galaxies into red galaxies is lower than 
the replenishment rate of blue galaxies, so that there is a 
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net decrease of red galaxies relative to the blue galaxies. 

However, the change in f re d with redshift does not give 
us unambiguous information on the actual changes in 
these rates. If we assume that the rate of the build-up 
of the red sequence is constant over the redshift range 
of ~ 0.8 to 0.2, we would then conclude that the rate of 
replenishment of blue galaxies (into the luminosity lim- 
ited sample) is decreasing from z ~ 0.8 to 0.6, and then 
becomes stable, or decreasing at a much lower rate, at 
z £0.6. Alternatively, if we assume that the replenish- 
ment of blue galaxies is constant, then the rate of trans- 
formation of blue cloud galaxies into the red sequence 
is increasing from z ~ 0.8 to 0.6, and become stable at 
z <0.6. The detailed picture is certain to be more com- 
plex. In a more general picture, the environment will be 
a major factor that affects these rates of flows. An im- 
portant effect to consider is the infall of galaxies into a 
high density region such as the parent halo of a galaxy 
cluster or group. We will discuss this effect in ^5.4[ in 
conjunction with the evolution of f re d in clusters. 

5.3.3. The fred of the Cluster Neighborhood Sample 

We have also derived f re d for the RCS-WSF cluster 
neighborhood sample for the z ~ 0.3 and z ~ 0.4 bins, 
which are plotted as open circles in Fig. [151 They show a 
clear difference due to the environments of the galaxies; 
the cluster neighborhood samples have f re d of ~ 0.8, in- 
dicative of high galaxy density regions. While we do not 
have the redshift range to examine the evolution of f re d in 
clusters, studies using similar techniques of large samples 
of clusters show a continuous decrease of f re d with red- 
shift, i.e., the Butcher-Oemler effect (Butcher & Oemler, 
1984). For instances, Loh et al. (2008), using a sample of 
approximately 1000 clusters from RCS1 at 0.4 < z < 0.9, 
show a steady decrease in f re A of about 0.4; while the 
spectroscopic sample of lEllingson et al.l (|2001[ ) indicates 
a steady drop of f re d from 0.9 to 0.7 at z=0.2-0.5. 

Comparing these results with the WiggleZ neighbor 
sample, there appears to be a rather different behavior 
in the change in f re d as a function of redshift for galax- 
ies in regions around star-forming galaxies from those in 
high galaxy density regions around massive halos. In- 
stead of a steady decline in f re d with increasing redshift, 
the f re d for the WiggleZ neighbors have a basically flat 
dependence on z, and show a significant drop only at 
z^0.65. Similarly, Iovino et al. (2010) show that galax- 
ies classified as being in a group environment have a 
much steeper dependence of their blue-galaxy fraction 
on redshift than those considered to be in isolated en- 
vironments. This difference in the evolution of f re d can 
be considered as a clear demonstration of the effect of 
environment on the evolution timescale of galaxy popu- 
lations. 

5.4. Galaxy Evolution and Environment 

The topic of galaxy evolution has been studied in 
galaxy clusters over many decades. With modern large 
surveys, the focus has extended to field galaxies. Here, 
we would like to use the term 'field' to refer to all galax- 
ies in a general blind field survey, such as SDSS or 2dF, 
where galaxy density ranges from that of isolated galax- 
ies to clusters. Studies using such 'field' samples have 
revealed a strong correlation between many properties of 



galaxies and environment (usually parametrized by lo- 
cal galaxy density). One example is that the fraction 
of star-forming (or passively evolving gala xies) changes 
strongly with local galaxy density ( e.g., iBaldrv et al.l 
[20061 : iLi et aTll200l llovino et al.ll2010l ). Other works fo- 
cus on star- forming galaxies. These star- forming galaxies 
naturally reside in regions beyond galaxy clusters. Their 
star-formation rates and colors are found not to have 
a str ong depende nce on environmen t (e.g . , iCarter et al.l 
l2Q01tlRines et al.l2QQ5l : fBalogh et a l. 2009; Cassata et all 
2007). Such conclusions appear to be different from the 
results from using all 'field' galaxies, and this is likely a 
reflection of the sample properties. 

Our WiggleZ galaxies are members of these star- 
forming galaxies, selected by their detectable UV flux. 
We use the term 'WiggleZ neighborhood' to refer to lo- 
cal regions around the WiggleZ galaxies to distinguish 
it from the general 'field' environment. The WiggleZ 
neighborhood is likely low-density environment regions 
in the large-scale structure. The galaxy sample gener- 
ated by counting excess galaxies in the WiggleZ galaxy 
neighborhood represents a census of galaxies in regions 
around star-forming galaxies covering a significant red- 
shift range. 

We also examine f re d in dense environments using the 
neighbors of the RCS-WSF sample of cluster galaxies. 
The much higher f re d value (~ 0.8 vs ~ 0.3) measured 
is an indication of the more rapid build-up of the red 
sequence in dense environments. Furthermore, it appears 
that the evolution of f re d for these two samples is also 
very different, with the WiggleZ neighbors having very 
moderate or no increase in f re d since z ~ 0.6. 

In the context of the discussion of Figure [16] in £15.3.21 
the flow paths and rates of blue and red galaxies on the 
CMD in cluster environments would be very different 
from those in the WiggleZ neighborhoods. Here, the 
major source of replenishment of blue galaxies in mas- 
sive dense halos is likely to be the infall of galaxies along 
the large-scale structure. These galaxies are then trans- 
formed into red galaxies within some time scale due to 
processes such as tidal striping, ram pressure, galaxy ha- 
rassment and interactions, which are generally associated 
with the quenching of star formation in galaxy groups 
and clusters. The continuous increase in f re d can also be 
in part attributed to the cosmological d ecrease in infall 
rate w ith time in a low Q m universe (e.g., Ellin gson et al.l 
I2Q03). 

There is an additional difference between these two 
samples in the shape of the GLF of the red galaxy pop- 
ulation. The Schechter function fit indicates that the 
red galaxies in dense environments have a GLF with 
a steeper faint-end. This can also be a result of the 
red sequence in cluster environments being in a more 
advanced build-up stage, or, the build-up of the faint- 
end of the red sequence in dense environment is more 
rapid than that in regions around star-forming galax- 
ies. The build-up of the faint-end of the red sequence 
in clusters has been traced by a number of studies (e.g. , 
Bell et al. 2004; Tanaka et al. 2005; Willmer et al. 2006; 
Stott et all 120071 : iDe Lucia et all 12009( 1. IGilbank et all 
(2008)) show some preliminary evidence that this build- 
up is faster in rich clusters than in poor clusters, which 
is consistent with the comparison here between clusters 
and low galaxy density regions. Another apparent differ- 
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ence in the GLF of the red galaxies in the two sample is 
that the red galaxies in the cluster sample show a signifi- 
cant excess to the Schechter function fit at bright magni- 
tudes, suggesting that there are relatively more massive 
red galaxies in dense environments, either due to initial 
galaxy formation history or an increased rate of mergers 
in the evolution of these galaxies. 

The possible difference in the evolution of M* of the 
red galaxy GLFs described in §5.2.2, while not of high 
statistical significance, also fits in with the scenario of the 
dependence of the evolution of the red galaxy population 
with environment. The lower value of Q (i.e., slower evo- 
lution of M*) for the Faber et al. (2007) 'Red' galaxy 
sample, which are in an environment denser than those 
from the WiggleZ sample, suggests that red galaxies in 
denser environments are likely older. This would be ex- 
pected if galaxies in dense environments such as clusters 
and group turn red and dead earlier in the history of the 
universe than those in low-density region. 

In summary, in the general scenario of galaxy evolution 
and its connection to the environment the WiggleZ neigh- 
borhoods can be considered as regions where there are 
minimal major environmental events affecting the evolu- 
tion of galaxies. These regions are likely the low galaxy 
density parts of the large cosmic structure where star 
formation is still occurring well into the current epoch. 
Here, galaxies are unlikely to be affected by environmen- 
tal influences that are associated with infall into a mas- 
sive halo such as a substantial galaxy group or cluster. 
The galaxies in these neighborhoods can be considered 
as primarily following a secular galaxy evolution path, 
controlled by their own nature at birth, with environ- 
mental effects playing a role over a much longer time 
scale. In these environments, only a small fraction (~20- 
30%) of the galaxies have turned red over the redshift 
range of 0.7 to 0.3. In a cluster region, where there is a 
continuous infall of galaxies from the lower-density envi- 
ronments replenishing the blue cloud, the effects of the 
large dark matter parent halo on galaxies would produce 
a very different mix of galaxy populations. In these re- 
gions, environment plays a dominant role, accelerating 
the quenching of star formation, and transforming the 
infallen galaxies into the red sequence over a relatively 
short time scale, producing a red sequence in a more ad- 
vanced evolutionary stage. 

6. SUMMARY 

We have probed galaxy evolution at 0.25 < z < 0.75 
using optical data from the RCS2 around ~41,000 star- 
forming spectroscopic galaxies from the WiggleZ project. 
Because of the complicated selection criteria in the Wig- 
gleZ survey, galaxies in the spectroscopic sample have 
discrete characteristics as a function of redshift. We 
therefore examine optical properties of galaxies within 
0.25 Mpc to WiggleZ galaxies using stacked color-color- 
magnitude cubes. The idea is to use the WiggleZ galaxies 
as markers, and assume that they and the surrounding 
neighbors are at the same redshift. By applying back- 
ground subtraction and stacking the net excess counts 
over a large number of markers, we are able to study the 
properties of the neighbors around the WiggleZ galaxies. 

We also examine how the optical colors of the neigh- 
bors correlate with the WiggleZ [OIIJA3727A equivalent 
width, NUV flux, and AGN activity indicators at differ- 



ent redshifts (Figures [7J [5J and [T0|). We find in general 
that the neighbor galaxies populate the same color-color 
spaces without significant dependence on these proper- 
ties of the markers, suggesting that the properties of the 
neighbor galaxies are not strongly affected by, or cor- 
relate with, the characteristics of the WiggleZ galaxies 
themselves, Thus, they can be used to study the evo- 
lution of the photometric properties of galaxies in low- 
density, star-forming regions over the redshift range of 
0.25 to 0.75. 

Our major findings are: 

(1) The majority of WiggleZ neighbors are blue galaxies 
which have a steeper faint- end slope and a faster evolu- 
tion term in 0,5 M*, than the red galaxies. 

The CMD of the WiggleZ neighbors shows the charac- 
teristic bimodal distribution of a red sequence and a blue 
cloud, with the latter dominating the galaxy population, 
containing 65 to 85% of the galaxies (depending on 
the depth of sampling and redshift). The GLFs of the 
two populations can be fitted with single Schechter 
functions, with the blue galaxies having a much steeper 
(a ~ —1.3) faint-end than the red galaxies (a ~ —0.4). 
Based on the three low-redshift bins (0.25 < z < 0.55), 
where the data are complete to 0,5 M r /=-19.0, we find 
no significant changes in the faint-end slope of the GLFs 
with redshift. There is significant evolution in M* for 
both the 'Blue' and 'Red' subsamples. We find that 
blue galaxies have a marginally more rapid evolution in 
M* over this redshift range with Q ~ —2.10, compared 
to Q ~ —1.59 for the red galaxies. While the Q values 
for both the red and blue galaxy samples appear to be 
steeper than the typical values in the liter ature, these 
evolut ion factors are similar to those from iFaber et al.l 
(|2007l ) when compared over the same redshift range. 

(2) The red galaxy fraction f re d in the WiggleZ neighbor- 
hood is approximately constant since z ~ 0.6, but drops 
atz>0.7. 

The WiggleZ neighbor galaxies have a red-galaxy 
fraction {fred) considerably smaller than that of the 
neighbors of the RCS-WSF sample of markers, at 
0.25 < z < 0.45, obtained the same way as the WiggleZ 
sample. The evolution of the WiggleZ neighbor f re d with 
redshift is modest; f re d can be described as basically 
fiat as a function of redshift with only a very moderate 
decrease up to z ~ 0.6. The average f re d value over 
the redshift range of 0.25 < z < 0.65 is ~ 0.30, similar 
to field galaxies at z ~ 0. A large drop to f re d~ 0.20 
is seen for the z ~ 0.7 sample. This drop is confirmed 
by extending the measurement of f re d to z ~ 0.8 using 
a sample with a slightly brighter absolute magnitude 
limit. Furthermore, this drop does not seem to be 
associated with the larger average luminosity of the 
markers at the higher z bin. The change in f re d with 
redshift in the WiggleZ neighborhood can be seen as 
either a higher rate of relatively bright (M^M* + 1) 
star-forming galaxies entering the luminosity limited 
sample at z ~0.7, or a decrease in the quenching rate of 
star formation at this redshift. 

(3) The comparison between the WiggleZ and RCS-WSF 
neighbor samples shows an environmental influence on 
galaxy properties and evolution, with the red sequence in 



22 



LH. Li et al. 



cluster environment being in a more advanced build-up 
stage. 

We examine the effects of environment on the galaxy 
population properties by comparing the WiggleZ neigh- 
borhood galaxies to that of the RCS-WSF neighbor sam- 
ple. Besides the expected and obvious difference in 
f re d values of the two samples (with f re d, at ~ 0.8, 
being much larger in the RCS-WSF sample), we also 
find significant difference in the GLFs of the red galax- 
ies of the two samples. The faint-end slope a for the 
red GLF of the RCS-WSF sample is considerably steeper 
(~ —0.7 vs ~ —0.4). This can be taken as the build-up 
of the faint-end of the red sequence being in a more ad- 
vanced stage in rich environments (such galaxy clusters 
and groups) than that in lower galaxy density regions 
around star- forming galaxies. Furthermore, there is also 
evidence that there are excess luminous red galaxies in 
the RCS-WSF sample. These findings point to the im- 
portance of environment in affecting the history of star 
formation in galaxies. Galaxies in cluster/group envi- 
ronments likely have suffered significant environmental 
events that rapidly shut down their star formation, turn- 
ing the galaxy red. Whereas in regions where star forma- 
tion is still prevalent, environmental events likely occur 
much less frequently and their effects spread over a longer 
time scale, delaying the build-up of the red sequence. 
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